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Abstract 


2 

The  paper  analyzes  the  spaces  55?(Q)  and  the  associated 
trace  spaces  on  the  boundary  dQ .  These  spaces  are  essential  in 
the  theory  of  the  h-p  version  of  the  finite  element  method.  The 
h-p  version  for  the  problem  with  nonhomogeneous  essential  and 
natural  boundary  conditions  is  analyzed.  Numerical  experimenta¬ 
tion  is  presented. 


Introduction 


1 . 

-  i 

There  are  three  versions  of  the  finite  element  method:  the 
h-version,  the  p-version  and  the  h-p  version.  The  h-version  is 
the  standard  one,  where  the  degree  p  of  the  elements  is  fixed, 
usually  on  low  levels,  typically  p  =  1,2  and  the  accuracy  is 
achieved  by  properly  refining  the  mesh.  The  p-version,  in  con¬ 
trast  fixes  the  mesh  and  achieves  the  accuracy  by  increasing  the 
degree  p  of  the  elements  uniformly  or  selectively.  The  h-p 
version  is  a  combination  of  both. 

The  standard  h-version  has  been  thoroughly  investigated  theo¬ 
retically  and  computationally.  The  literature  here  is  overwhelm¬ 
ing.  To  date  there  are  over  two  hundred  monographs  and  conference 
proceedings  [18]  and  new  monographs  and  proceedings  are  continuous¬ 
ly  appearing.  There  are  many  programs  of  research  and  commercial 
type  available  (e.g.  see  [18]). 

The  p  and  h-p  version  is  a  new  development  and  it  is  very 
successfully  used  for  solving  elliptic  equations,  especially  in 
the  field  of  computational  mechanics.  The  first  theoretical 
results  were  published  in  1981  (see  [2], [10]).  There  is  only  one 
commercial  code  based  on  the  p  and  h-p  version  of  the  finite 
element,  the  program  PROBE  of  NOETIC  Technologies  (St.  Louis,  MO). 

PROBE  deals  with  two-dimensional  elasticity,  stationary  heat 
problems  and  thermoelasticity  problems.  The  code  for  the  three- 

dimensional  problems  will  be  released  in  1988.  PROBE  presently  is 
utilizing  Is  p  s  8.  There  is  also  commercial  code  FIESTA  for 
solving  three-dimensional  elasticity  problems  using  1  <  p  *  4.  A 


research  code  STRIPE  developed  by  Aeronautical  Research  Institute 


of  Sweden  has  the  p  and  h-p  version  features  for  three- 
dimensional  problems  and  is  using  2  <  p  <  12. 

For  the  survey  of  the  today's  state  of  the  art  and  recent 
progress  we  refer  to  [ 1 ] , [ 2 ] , [ 8 ] , [ 14 ] , [ 19 ]  where  also  additional 
references  can  be  found. 

The  success  of  the  h-p  version  is,  among  others,  based  on 
the  fact  that  the  elliptic  problems  of  the  structural  mechanics 
are  usually  characterized  by  piecewise  analytic  data  (boundary, 
coefficients,  boundary  conditions).  This  implies  then  that  the 
exact  solution  is  analytic  (or  piecewise  analytic)  with  singular 
behavior  of  precise  character  in  the  a-priori  known  areas  as  for 
example  in  the  neighborhood  of  the  corners  of  the  domain.  We  have 
shown  in  [4], [5]  that  this  class  of  solutions  can  be  very  accu¬ 
rately  described  in  the  frame  of  countably  normed  spaces.  We  have 

2 

denoted  this  space  by  8^(0).  If  the  solution  belongs  to  the  this 
class  then  we  have  shown  in  [6], [13]  that  the  finite  element  solu¬ 
tion  converges  exponentially. 

The  present  paper  elaborates  on  the  characterization  of  trace 

2 

spaces  of  the  function  u  e  8^(0)  and  9ives  precisely  verifiable 

necessary  and  sufficient  conditions  for  the  input  data  (Dirichlet 

and  Neumann,  conditions,  right  hand  side)  which  guarantee  that  the 

2 

solution  belongs  to  8^(0).  In  the  previous  paper  we  did  address 
the  h-p  version  for  the  problems  where  the  essential  (Dirichlet) 
boundary  conditions  could  be  satisfied  exactly  by  the  finite  ele¬ 
ment  solution.  In  the  present  paper  we  design  and  analyze  the  way 
how  to  deal  with  nonhomogeneous  essential  boundary  conditions  in 
the  full  generality.  We  show  that  the  performance  of  the  method 


is  the  same  for  general  essential  conditions  as  for  the  natural 
ones.  In  section  2  we  give  the  preliminaries  and  basic  defini¬ 
tions.  Section  3  defines  the  model  problem  of  second  order  ellip¬ 
tic  partial  differential  equations.  Section  4  introduces  the 

2 

spaces  of  traces  of  u  e  $^(Q)  on  r.  It  shows  also  that  the 

2 

function  in  the  trace  spaces  can  be  extended  into  This 

section  gives  some  of  the  major  results  of  the  paper.  Section  5 
defines  the  finite  element  method,  its  h-p  version,  characterizes 
the  meshes  and  elements  under  consideration  and  defines  how  to  deal 
with  nonhomogeneous  boundary  conditions.  Section  6  is  analyzing 
the  convergence  of  the  method  and  proves  that  the  rate  of  conver¬ 
gence  is  exponential.  Finally,  Section  7  brings  numerical  exam¬ 
ples  which  show  that  the  theoretical  results  having  an  asymptotic 
character  are  applicable  in  the  wide  range  of  practical  accuracy. 


is  link! 


domain  Q  a  straight  polygon.  Otherwise  we  will  speak  about  a 


curvilinear  polygon.  If  0  <  <  2 n,  i  =  1,...,M,  we  will  speak 

about  a  Lipschitzian  domain.  Let  us  assume  that  r  = 

where  r(0)  =  U  f  ,  r(1)  =  r-r(0),  f(1)  =  U  f.,  where  Q  is 

ieQ  1  ieQ' 

some  subset  of  the  set  <1,2,...,M}  =  /ft  and  Q'  =  ft-Q. 

We  have  assumed  for  simplicity  that  Q  is  a  simply  connected 
domain.  The  results  we  are  presenting  here  are  also  valid  when  0 
is  n-connected,  bounded  domain  and  its  boundary  is  composed  by 
n-curves . 

Denote  I  =  {x|-l  <  x  <  1},  we  also  will  write  I  =  {x  ,x  | 
-1  <  Xj^  <  1,  x2  =  0}  c  ip  when  no  misunderstanding  could  occur. 

By  L9(Q),  L  (Q),  L_ { I ) ,  L  (I)  the  usual  spaces  of  p- 
integrable,  1  <  p  <  ®>,  functions  on  Q  or  I  are  denoted.  By 
Hra(Q),  Hm(I),  m  2  0  an  integer  we  denote  the  usual  Sobolev  space 
of  functions  with  square  integrable  derivatives  of  order  <  m  on 
Q  (respectively  I).  The  space  Hm(Q)  is  furnished  with  the 
usual  norm 


|U|! 


Hm(Q  ) 


I 


!!  Daul!  2 
L  - 


(Q) 


0<  |a | <  m 


where  a  *  (ai'a2^'  ai  '  0  integer,  i  =  1,2,  |a| 


=  a1+°2  and 


_a 
D  u 


d  u 


axj'ax®* 


=  u 


O  t  ft  t 
X1  x2 


Further  we  let 


u|  =  II  I  Dmu  |  ||  ...  , 

Hm(0)  L2(Q) 


ID 


m 


u|2  =  |Dau|2. 


I  a  I  =m 


v*  .■>  , 


As  usual  we  shall  write  H  (0  )  =  L 2(0), 

hJ(Q)  =  {u  e  H1 (0 } | u  =  0  on  r(0)}. 

In  the  analogous  way  we  define  Hm(I)  with  Dku  = 

dxK 

By  rj(x)  =  dist(x.Aj)  =  Ix-A^l,  x  e  0,  j  e  M  we  shall 
denote  the  Euclidean  distance  between  the  point  x  and  the  vertex 
Aj,  r1(x)  =  lx+1^  r2(x)  ^Ix-l^  x  e  I.  Let  3  =  (3  . . . / #  M )  (res¬ 
pectively  3  =  (3^,3^))  be  an  M-tuple  of  real  numbers  0  <  3^  < 
1,  i=l,...,M.  We  will  write  a ^  <3  <  (respectively  3  <  3) 

if  ai  <  <  a2  (respectively  3^  <  3^),  i  =  For  any 

k  integer  we  shall  write  3+  k  =  {3^+k,  .  .  .  ,#M+k}  and 


*tf+k(x)  =  n  I  rj,(x) \3i  +k,  xe  Q 


*s+1_(x>  =  n  lr,(x)  |/?i+k,  x  e  I. 
i=l  1 

By  C^(Q),  c^(Q),  C^(I),  C^(I),  j  >  0  integer  we  will  denote  the 

set  of  all  functions  with  continuous  j-derivatives  on  0,0,  I, 

I,  furnished  with  the  usual  norm  :!  •  I!  .  ,!!•!!.  .  Let 

CJ(0)  CJ(I) 

in  ^ 

'  (Q )  ,  m  2  t  >  0  integers  be  the  completion  of  the  set  of  all 
infinitely  differentiable  functions  under  the  norm 


bJ-'io) 


,Ull^-l(01  *  I  "W<ID'e"“'»L2<0) 

H  101  lai-e  2 


for  f  >  1 , 


=  V  !!♦„  .,JDau!  !!?  _ 


If  m  -  £  =  0  we  shall  write  '  =  L^(0).  Analogously  as 

before  we  define 


4'^Q)  |a|-€ 


L2(0)* 


In  the  similar  way  H™'  (I)  is  defined 


!!u!!  2_  ,  =  !!  u!!  2 e  . 

Hj'  (I)  n  x(D 


f°r  f  ! 


lim  '  0  /  T  \ 

«  ( 1 ]  k=0 


-  X!w°0“'bl 


"l2(d  • 


Further  we  introduce  the  space  8^(0),  l  >  0  integer  which  will 
play  an  important  role  in  this  paper: 


*8  <°  I 


<ulu«  any  k  ,  t,  Daul  E  ^  (Q  , 

s  Cdk-^(k-€)  !  ,  |a|  =  k,  C  >  0,  d  >  1 


independent  of  k} . 

If  we  wish  to  underline  the  dependence  on  d  we  will  write 

t 

Analogously  for  f.  >  0  integer 


8^(1)  =  { u  I  u  e  hJ'€(I),  any  k  >  (,  lit 


,,  (k)., 

fi+k -lu  "  l2(I) 


Cd  (k-£)!,  C  >  0,  d  >  1  independent  of  k> . 


Further  for  j  =  1,2, 


«rj{0)  -  {u  e  H^'J(0)  ||Dau|  (x)  <  Cdkk!  |*k+/J_j  +  1(x)  r1 

|a|  =  k  =  j-l,j,...,C  >  0,  d  >  1  independent 
of  k)  , 


<r^(I)  =  iue  H^J(I)||u'"'|  <  C|4k+/?-j  +  l/2(X)  1  *d"k! 


k  >  >  0,  d  >  1  independent  of  k} . 


_  1/2 

Let  y  e  U  r,.  Then  we  define  H  (y)  (respectively 

iepc^i 


Hk  ^  ^  (y )  ,  k  ^  €)  ,  k  ^  €  >  1  integers  as  follows:  for  any 


<*>  e  Hk  1//2(z)  (respectively  Hk  1//2(^))  there  exists  f 


Hk(Q),  (respectively  Hk'^(Q))  such  that  f  j  ^  =  <p .  We  define 


then 


!Hk-l/2,r,  ( respectively  S*S  k.1/2<<.1/s 

«  W  J  1 


=  inf  ||f ||  .  (respectively  iifll  .  f  ). 

f|  r^p  hk(0)  hx'c(Q) 


H?  (Q) 


By  » 


€-1/2 


(/),€>!,  we  will  denote  the  set  of  the  traces  on 


of  functions  from  the  space  »  (ft ) . 


Let  rd  be  an  edge  of  Q,  then  by  the  assumption  there 
exists  a  one  to  one  mapping  nu  of  I  onto  which  is  analy¬ 


tic.  If  Fq  is  a  straight  line  then  we  shall  assume  that  i: 


the  linear  mapping.  Let  u  be  defined  on  U(x)  =  u(mi(x) 


be  defined  on  I.  Then  we  define 


Hm(ri)  =  (u | U  e  Hm(I) } 


Hull  =  HUH 

Hm(ri)  Hm(i) 


m  l  £  f 

In  the  same  way  we  define  the  spaces  H™'  (f^),  e^'(r^) 


Let  us  remark  that,  as  we  defined  it,  :!  •  I!  „  depends  on  the 

Hm(F.) 


mapping  mA ,  i.e.,  it  depends  on  the  parameterization  of  the  arc 


r i .  Nevertheless  the  space  H™'  (f^)  does  not  as  well  as  ff^(f  ) 


(see  Lemma  4.6)  but  35  ^  (f  ^ )  could  be  dependent  on  m^ .  Let  us 


state  now  some  lemmas  which  will  be  used  later. 


Lemma  2.1.  We  have 


H *'2(Q)  c  C°(Q) 


with  the  continuous  injection. 
See  Lemma  7  of  [3]. 


Lemma  2.2.  Let  u  e  '  (Q ) .  Then 


(2.1) 


ii  I  D  u!  §  l|  <  Ciluii 

3-l]L2(Q )  "  "h2.2(0) 


(ii)  Let  u(Aa)  =0,  i  =  1, . . . ,M.  Then 


(2.2) 


IUV2"l2(Q)  S  C!IU"h2,2(Q) 


See  Lemma  8  of  [2].  ; 

Lemma  2.3.  ®^(«)  c  ^(Q)  and  **(Q)  c  (« )  ,  0  </?+£<  1,  f  > 

0  arbitrary.  : 


See  Theorem  2.2  and  2.3  of  [6]. 


Lemma  2.4.  Let  u  e  8^(0),  j  >  o,  then  u  is  analytic  on  0  - 


U  A,  . 

1  =  1  1 

Lemma  2.5.  Let  r  0  1  and  F(x) ,  0  <  x  <  00  is  defined  by 


F(x)  =  j  ^{t)dt  for  r  >  1 


F(x)  =  f(t)dt  for  r  <  1. 


Then 


3.  The  model  problem  and  its  properties 

Let  0  be  the  curvilinear  or  straight  polygon  and  L  be  a 
strongly  elliptic  operator 


C.  £ 

U)  *  ■  X  (ai.3(x,ux ,»*.  *  Zbilx)ux. 

■:  _  i  X  J  2  —  *  X 


+  c  (x)u 


i.  j=l 


where  a.  .  (x)  =  a.  .  (x),  b. (x) ,  c(x)  are  analytic  functions  on 

X  *  J  J  0  x  1 

Q  and  for  any  ^i'^2  e  ^  and  any  x  e  Q  let 


X  •i.j'i'j  *  V'X’ 


i.  3-1 


with  vQ  >  0. 


Let  B(u,v)  be  continuous  bilinear  form  on  Hq(Q)*Hq(Q) 


r  2 

iu'v|  =  X  a1.3Ux,vx,  *  X  biux, 

JQ  ''i ,  1=1  J  i=l 


v  +  cuv  dx. 


We  assume  that 


inf  sup  |B(u,v)|  >  (i  >  0 

lull  =1  ||  v  !i  =1  1 

H  («)  H  (Q) 

ueHj(Q)  veHj(Q) 


and  for  any  v  e  H q(Q),  v  x  0 


sup  I B(u, v) |  >  0 . 
lull  ,  =1 

H1(0) 

U^HQ(0) 


Asstime  now  that  g[€]  e  *2/2_€<r  (€)  )  ,  t  =  0 , 1 ,  f  e  ?J(0)  and  con¬ 
sider  the  boundary  value  problem 


(3.1a) 


Lu  =  f  on  Q 


(3.1b) 


L  u  J 

u  =  gl 


__  rlW 

on  I 


(3.1c) 


du 

dn 


-  gm 


on 


1) 


where  we  denoted  by  n  the  conormal  in  the  usual  sense.  The 

c 

solution  of  our  problem  is  understood  in  the  usual  sense.  Then  we 
have 

Theorem  3.1.  There  exists  unique  solution  uQ  e  H1 (Q )  of  the 
problem  (3.1).  See  Lemma  3.1  of  (4].  : 

Let  us  mention  some  theorems  addressing  regularity  of  the 
solution  uQ. 

Theorem  3.2.  There  exists  0  s  tf.  <  1,  i  =  1, . . . ,M  depending  in 
the  problem  (i.e.,  operator  L,  <o.,  etc.),  such  that  if  f 
«°(Q),  €  »|/2_'f(r  )  ,  t  ■  0,1,  B  <  S  <  1  then  uQ  e  **(0)). 

Proof  is  given  in  [3]. 

Theorem  3.3.  Let  Q  be  a  (curvilinear)  polygon  (instead  of 
straight  polygon  as  in  Theorem  3.2)  and  let  then  assumptions  of 
Theorem  3.2.  hold.  Then  uQ  e  ff2(Q). 

Proof  of  the  theorem  is  given  in  [6], 

We  have  seen  in  [6],  [13]  (see  also  sections  5  and  6)  that 

when  the  solution  u  of  the  problem  3.1  belongs  to  the  class 
2 

'Pj(O)  then  the  h-p  version  of  the  finite  element  method  con¬ 
verges  exponentially. 

Theorems  3.1  and  3.2  show  that  it  is  important  to  develop 

3/2— C 

practical  characterizations  of  spaces  *  (f  )  ,  *’  =  0,1.  which 

can  be  easily  used  in  concrete  cases  to  verify  whether  the  input 

m 

data,  i.e.,  g  belong  to  the  desired  space.  We  will  elaborate 
on  it  in  the  next  section. 

12 


4  . 


Traces  and  extensions  of  weighted  Sobolev  spaces.  Characteri¬ 
zation  of  the  spaces  8^2  ^(F  ) 

In  this  section  we  will  elaborate  on  the  characterization  of 
the  space  8^  (T  )  ,  l  =  0,1  which  leads  to  an  easy  verification 
in  the  concrete  cases  of  applications. 

Lemma  4.1.  Let  13  =  (!3 1,fi  2)  ,  0  <  fi  <  1/2  and  g  €  H  *  ' 1  ( I )  .  Then 

(i)  g  e  C°(I)  and  :igll  _  <  Ci!g!! 

c  (I)  H  J  '  <  I ) 

(ii)  I  g(x)  -g(  -i  )  |  ,  ci1/M(x)llgl| 

"/?  ' 1 ' 

1  g(x)  -g(  i  )  |  ,  «1/M(x)Hgil 

“13  ( X) 

where  C  is  a  constant  independent  of  g(x)  (but  depends  on  3) . 
Proof .  Obviously 


( I ) 


|g(x)-g(t)  I  s  I  g'(T  )dr  | 

t 


(4.1) 


r  x  , 

g'2(T  )**(T)dT 

1/2 

r  f  X  -| 

(♦,,(*)  r2dr 

-  t 

-  t 

1  1.1  (♦^(r))'2dr 

Hi*  *i)Ut 


which  shows  that  g  is  continuous  on  I.  Using  the  imbedding 
theorem  on  (-1/2, 1/2)  =  I'  we  have 


(4.2) 


I  g  (  0 )  |  s  Cllgli  s  Cllg'! 

HA  ( I' )  '  1  ( I ) 


and  we  get  immediately 


Further  (4.1)  immediately  leads  to  (ii). 


Lemma  4.2.  Let  /?  =  (3^,3^),  1/2  <  ^  <  1  and  9 


2  2 

H  3  '  (  I ) 


Then 


(i)  g  e  CU(I)  and  ilgll  .  <  C !| g||  .  _ 

r*u  /  t  \  '  * 


CW(I) 


«  ( I ) 


(ii)  lg(x)-g(-l)  I  <  C#3/2_/?  (x)  ilgll  2  2 

H3  '  ( 1 ) 


|g(x)  -g(  1 )  |  <  C#3/2_/?  (x)  Ilgll  2  2 


H/3  '  I ) 


where  C  is  a  constant  independent  of  g(x). 


Proof .  Using  (4.1)  we  get 


I  g(x)  -g(  t )  I  s  |  g'(r)dr| 


^  g'2(T)4"^dr 


5  l!g'*l-/J!!L0(I)  *I-«<T>dT 


s  Hg'-a'ioni^s^d, 


+  I  g'  ( o )  |  ii  * 


1-/3  1  L2  (  I ) 


5  C[|g'(0)|  ♦  (I)] 


S  Cllgll  2<2 


h;*‘(I) 


In  the  last  inequality  we  used  Lemma  2.5  and  the  fact  that  1/2  < 


3  <  1 .  The  lemma  now  follows  Immediately. 


Lemma  4.3.  Let  g  e  d(I),  0  <3  <  1.  Then  for  k  >  1 


'  ' ‘  '  k-l/2+ j  '  '  1 ' 

where  y >  1  is  independent  of  g,k,d,  and  C  depends  on 
but  is  independent  of  g,k. 


Proof .  Let  I'  =  (-1/2, 1/2).  Then  for  any  k  i  1  we  have 

"g(k)!i  ,  <  C(f  (1/2)  )~k_*k!dk 


H1  ( I' ) 


where  3  =  max(  $  ^  ,  <3^  )  .  Hence  by  the  imbedding  theorem 


|g(k)(0)l  <  Cdkk ! 


where  dJ  2  yd,  y  >  $  (1/2)  >  1.  Further,  for  k  <  1 


g(k)(x)|  <  lg(k)(0)|  +  |  g(k+1)(t)dt! 


5  lg{k) (0) I  + 


(g(k+1) (t) )2*^+k(t 


<  C(d, )*k! (♦ 


k-1/2+3 


4.4. 

Let  g  e  8  * ( I)  ,  0  < 

(  <  1 . 

Then 

9  -  *  5  < 1 1 

4.5. 

Let  g  e  8*(i)  ,  0  < 

t  <  1 . 

Then 

for  k  ■  2 

( k) 

ig  (x)  1  <-  c(*k.3/2+<(x) > 

- 1  ,k. 

djk ! 

and  g  e  <?,..(  I )  . 


*"»Tr^r 


.1  1.1^1 


Lemma  4.6.  Let  K  -  m(x)  be  a  one  to  one  map  of  I  onto  I, 


m(x)  be  analytic  on  I  and  !m'(x)|  >  0,  x  e  J.  Assume  that 


g  e  e~(I),  j  =  i,2,  and  define  v(x)  =  g(m(x)).  Then  v  e  ffj(l), 


j  »  1,2 


Proof .  Because  m(x)  is  analytic  on  I  it  can  be  extended  into 


the  complex  plane  €  on  I§  =  {z  *  x+iy|-l-6  <  x  <  1+6, |y|  <  6), 


6  >  0,  m(z)  is  a  one  to  one  mapping  of  onto  o  i  6'  > 


0  and  |m' (z)|  >  an  >  0,  z  e  I .  .  Let  now  j  *  1  and  x. 


Then  for  k  >  1 


IOtk,<V'  !  Cl*k-l/2«'X0>>"ld^! 


and  the  series 


g'(x)  -  £ 


«'**"lx0|lMo|kJ! 


♦(x0) 

is  absolutely  convergent  for  |x-xQl  s  a — g — ,  a  <  1,  Hence  also 


g'(z)  *  ^g(k+1,(x0)(z-x0)kj! 


--1  ♦(x0> 

converges  and  |g'(z)l  *  C*/?  +  1/2<X0)  for  1  Z“X0 1  '  a — where 


C  Is  independent  of  xn .  Hence  g(z)  is  a  holomorphic  function 


and  v(z)  ■  g(m(z))  is  holomorphic,  too.  Using  Cauchy  theorem  we 


get  Immediately  that  for  k  2  1 


I v<lt|  1x1  I  !  Cd2ik-l/2«lx,k!- 


Obviously  v{x)  e  '  (I).  In  quite  a  similar  way  we  prove  the 


statement  for  j  =  2. 


Remark  4.1.  Lemma  4.6  shows  that  the  space  <T  ^  ( I )  is  invariant 
with  respect  to  an  analytic  mapping.  Using  the  formula  of  the 
n*" h  derivative  of  a  composite  function  (see  formula  0.430  of 
[15])  we  can  also  show  that  (I)  is  invariant  space  with 
respect  to  an  analytic  mapping  m(x)  as  in  Lemma  4.6. 

Let  T  be  an  analytic  arc.  Then  we  could  define  the  spaces 
and  $(^(r)  with  respect  to  the  length  instead  as  we  did  in 
section  2  using  a  specific  mapping.  These  two  definitions  are 
then  equivalent  by  lemma  4.6.  and  Remark  4.1. 

Lemma  4.7.  Let  M(x) ,  x  e  9^,  M(x)  =  (M  (xJ.M-fx))  is  a  one  to 

one  mapping  of  0  onto  0  and  I J  1 1  <  /j  on  0,  where  J  is 

the  Jacobian  of  the  mapping.  Assume  that  M(x)  can  be  anaiyti- 

2 

cally  extended  on  =  (x  e  !?  !dist(x,0)  s  6  }  so  that  it  is  one 

to  one  mapping  of  0&  onto  0,0  o  0.  Let  u  e  ff^(0),  j  =  1,2, 
v(M(x) )  =  u(x) .  Then  v  e  ff^(Q).  The  proof  is  quite  analogous 
as  of  the  Lemma  4.6  only  we  have  to  apply  the  theory  of  two  com¬ 
plex  variables. 

Lemma  4.6.  Let  ge  0  <  '<  <  1,  j  =  1,2.  Then 

g«e  *^(I),  0  <  'T  <  l, 

/3 

Q  =  n+f.  ,  c  >  o  arbitrary. 

Proof .  Let  us  consider  only  the  case  j  =  1 .  The  case  j  =  2  is 
analogous.  Because  for  k  >  l 

|g(k)(x)|  <  Cdkk!  (♦k+l?_1/2(x)  )-1 


we  get 


(g(k)  (x)  )  (x)dx 

*  _  1 


<  Cd2k(k!)2  i2 


( x)  dx 


ft -ft -1/2 


2k  2 

5  C(c)d  (k!)  . 


We  see  that  Lemma  2.3  has  a  completely  analogous  version  for  the 
relation  between  $2(I)  and  ®2(I). 

k+2  2 

Theorem  4.1.  Let  u  e  '  (Q )  ,  k  >  0  and  r  be  a  straight 

line  edge  of  0  and  u|r  =  g..  Then 

i 

(i)  For  1/2  <  <  1  and  k  i  0 


»i  *  !ri,i|ri''  k  * 


K. 


>  0,  fti  e  (/?i  +  j_1-l/2.1)  ,  j  =  1,2 


!'  q . !!  .  <  cdk!!  u!i  ,  _  _ 

i  Hk+1 »l(r  )  Hk+2, 2(0 ) 


i 


with  C  independent  of  k  and  d  >  1. 

(li)  For  0  <  «i,/?i+1  <  1/2,  k  *  l 


SI  e  H  (P  i  )  , 


,k+l ,  2 


gA  e  H?  ( T  ± )  ,  fti  j  e  (/3i  +  j_1  +  l/2,  1)  ,  j 


!!  a±  it  ,  <  Cilusi  2  2 

H  ( r  A  )  H*'*(Q) 


(iii)  If  u  ■=:  'B2(G)  and  1/2  <  +  1  <  1.  then  gA  -- 

**(ri)'  ?i,j  ~  Hi  +  J-l'1/2'1/2)  '  j  =  1,2‘  If  0  <  *  i '  *  i  +  1  <  1/2 
;J  i 

then  gj  «  *?  <r±)  .  ') i  }  e  (^^.^1/2,11. 

'  i 

Proof ■  Without  any  loss  of  generality  we  can  assume  that  !  <  =  !" 
and 


r :  =  {x1,x2ix1  -  i ,  x2  =  o},  A1  =  (-1,0),  a2  =  (i,o). 


aku 

Let  k  >  0  and  v,  =  — r-*.  , .  Then  for  k  2  2 

K  dXj  K 


v  i. 

k  2 
H  (0) 


C[|:  ( D2  1-g)* 


axk,lk+-3  L2(0) 


+  k  (0^  ^ _ H)  $ 

axk  k+*’1  l2(0) 


+  k2  (^V*  1 

axk  k+'J-2  l2(0)J 


<  ck  u ■; 


Hk+2 , 2(Q  j 


Using  Lemma  2.2  we  get  for  k  =  1 


Vl'  2 
H  (0) 


C  u 


H2,2(0) 


Because  of  Lemma  2.1  u  -  C°(0)  and  hence  vQ(A^)  =  °<  *  =  1.2 
Hence  using  Lemma  2 . 2  we  get 
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.* 

V 

V 
’/ 

V 
■y 
< 

,N 

p> 


',1 

;o 

\A 

m 

% 

K 

$ 


r- v.v;.- vj* 


<  c !!  u!! 


4.*  (.4  j.|  |.|  l.«  4,<  4.«  aj> 


I-I'I-I  hUiaiXLfcbJjl, 

li 


v0!  2 
u  m  (Q  ) 


2  2 

a*'  (Q) 


and  hence  for  all  k  *  0 


(4.3) 


k»2,„.  S  C«lt+1>  l!"!„k+2.2 


H  (0  ) 


where  C  is  independent  of  k.  Therefore  by  the  imbedding  theorem 

vk  e  C°(Q) ,  k  >  1. 

Let  us  show  now  that  v^ ( A^ )  =  0,  i  =  1,2,  k  *  1.  Assume  on  the 
contrary  that  vk(A^)  >  Then  because  e  C  (Q )  we  have 

vk(x)  >  c  >  0  for  Ix-AjJ  <6,6  >0. 


Hence  for  k  *  2 


r  „  ak  . 

.2  ,d  u.  2.„ 

00  >  *k+/?-2 (7Tc  dX 

Q  X1 


2  2 

>  e  *l2dX  *  °° 


*-2vkdx 


where 

=  Q  n  { x 1 {x-A1 |  <  6 } 

and  we  have  the  desired  contradiction.  For  k  =  1  we  use  Lemma 


2.2  and  get 


If  u  c  9ff(Q) 


*U  2  2  .  2  [  .2  , 

*  >  <5x7l  f  *-2dx 

J°  °s 


then  we  get  from  (4.3)  for  k  >  0 


lv.  li  2  *  CdV'k !  . 

*  H  (Q) 


WRWW5 


tol  I— • 


We  have  g^(x  )  =  ^  |  _  ,  k  i  0.  Then  g^(x  )  = 

1  3x^  1 

”1  “1  — 1 

*  k+tf  (x)vk(x)  |r  =  $k+,5  (x!  >vk(xl J  where  we  wrote  §k+:J(x1)  and 

v^fx^  instead  of  *^+-3  ^xi  •  °)  and  v^x^O).  Assume  first  that 

w  _  rt 

<  3  ,  >3  <  1 .  Let  dfl={  min  dist(A.,r  )r  .  Then  we  have 

j=3 , . . . M  J 

for  X  r  ,  $(  xJ  )  s  <t(x1)dC)1  and  hence  for  j  =  1 , 2 ,  .  .  .  ,  k+1 

*2  -  I g( ^ (x  ) I 2dx 

J  3-1+/3 
l 

,1 


2 

s  Cj  { 


-1 


-  *  cdo2Jj2< 


1  2  ~  ? 
[IV.  ,1 2tf 

-1  V’ 


,  1 2:2 

+  I  V  .  ,  I  # 

J  -1 +3-3 


JdV 


Using  Lemma  2.5  the  fact  that  j  =  l,...,k+l,  v^_1(Ai)  =0,  i 
1,2  and  that  fl-'^  +  l  >1/2  we  get  for  some  d1  <  1 


$ 


r 1  J-1+"i 


lg(j) (Xj) I 2dx:  <  Cd~ 2 ^ 


-l 


v'.  I  2$  ?  dx  . 
J-l  ?i-j  1 


By  (4.3)  and  the  imbedding  theorem  we  have  for  1  <  p  <  *>  and 
j  =  1 , . . . , k+1 


?Vj-ieX.  (I)  5  c(P)‘  vj-i:  „2|q| 


c  j  2 '  u 


H j  +  1 ’ 2 (Q  ) 


Hence  for  j  =  l,...,k+l,  because  >  -1/2  we  get 


1  1 

*2  -  !g( j,<x  ) |2dx  s  Cd~2k{  (#?  )pdx) 1/Pj v'  2 

J  1  1  J_!  J_1  L2qU 


,2  ,p,,l/p,  ,  2 


_,-2k  2  „  ,-2k  2 

-  ^di  o  5  ^d2  U  k+2  2 

1  J  1  H  (Q)  H*  '  (O) 


Because  by  Lemma  2 . 1 


'ul  0  -  5  C'U  2  2 
C  (0)  Hj'*(Q) 


we  get 


’  9  ■  r  /  r  \  4  C  u  2  o 

L2l  1}  H,'2(Q) 

iS 


Hence  we  have  proven  (i)  and  (iii)  for  1/2  <  <  1  and 


k  2  0.  Assume  now  that  0  <  <  1/2.  We  will  proceed  analo¬ 


gously  as  before.  For  j  2  2  we  have 


*2  .  |g(j) (x  ) |2dx 

r,  ^2+/3i 


*  Cj' 


2  12 


2.-2 


§  r 1 v'  it  +t  1  v 

_a  j-2^1  J”1  *+j-1  ^ J"1 


-]dx1 


<  Cd  2j  t2  .  ( v'_  (x  ) ) 2dx 

J_1  -1 +B1-3  J 


when  we  once  used  Lemma  2.5  and  the  fact  that  -1+'?^-?  >  -1/2 


Hence  using  (4.3)  and  realizing  that  -1+3 ^-3  >  -1/2  we  get  analo¬ 


gously  as  before  for  j  =  2,...,k+l 


r 1 


|g(  ^  )  (x  )  |  2dx  <  cd21C  u‘  2k+2  2 

1  1  ^  h2  '  (Q ) 


Let  us  prove  now  that 


m  •  »  ,  »  *  .  »  .  "  +  .  *  .  •  .  *  ,  ■  ««•_***  *  *  «  •  ,  «  «  1 


.1 


<  Cl  u  if  >  i 

1  °  u  I, +  2  2  '  K  -  1  ' 

H  (F.)  Hj  *'*(Q) 


We  have  VgfA^  =  vo  ^2  ^  =  0  and  hence 


2  -2 
g'  dx  <  CdQ 


-1 


[*..  2  I  vlt  2  +  (  v  1  V  2  ]dx  <  Cd  2 


-1 


3+1 


-1 


~ -2  2 
4  l  v'  I  dx 


where  we  have  once  more  used  Lemma  2.5.  Because  0  <  *  <1/2  and 


vo  vn  !  C<p|i'’°:h2(0 


<  C;  U 


2  2 
'  (0 


we  proceed  as  before  and  (ii)  and  (iii)  follow  easily. 


Remark  4.2.  It  was  essential  in  the  proof  of  Theorem  4.1  that 


>3i,j  ~  (/3i  +  j-i_1/2  '  1 )  respectively,  ^  e  (/3i  +  j-1  +  1/2  '  1  )  • 


i.e.,  of  the  open  interval.  The  proof  does  not  hold  for  the 


closed  interval.  It  was  assumed  in  Lemma  4.9  that  the  edge 


0 


of  the  domain  was  straight.  Let  us  assume  now  that  F  ^  =  m(I) 


where  m  =  (<p,v)  are  analytic  functions  on  I  as  given  in 
section  2.  Then  we  have 


Lemma  4.9.  Let  the  edge  of  the  domain  be  analytic.  Then  the 


part  (iii)  of  Theorem  4.1  holds 


Proof  .  By  Lemma  2.4,  u  e  5^(0).  Let  M(£  )  =  (<p  (t  )  ,  w  ( t  )  )  ,  t  -  I 


be  the  mapping  of  I  onto  F 1 .  Then  we  define 


M  {t  ,rt  )  =  <s>  (Y  )  -nw'  (Y  )  ,  M  (Y  ,n  )  =  i«  (Y  )  +ruo'  (t  )  . 
1  .....  2  " 


Then  the  mapping  M  (Z,ri)  =  ( M1  ( <  ,  r? )  ,  M2  ( <  ,  r? )  )  is  analytic  on 


Ic  =  (r.,T)\-l-f>  <  Y  <  i+ft,  |^|  <  6),  6  >  0,  |  J  |  <  o. 


-1 


<  <>  on 


(where  J  is  the  Jacobian  of  the  mapping)  and  maps  I,  onto 


the  (open)  neighborhood  S  of  F  ^ .  Denoting  0  =  0  S  , 
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k 


T  =  M  (Q  )  ,  we  see  that  v(x)  =  u(M  (x) )  is  defined  on  T, 

2  2 
and  v  e  <F  ^  ( T )  by  using  Lemma  4.7.  Hence  v  e  $^+f(T),  >.  >  0 

arbitrary,  by  Lemma  2.3.  Hence  for  1/2  <  <  1  we  get  by 

(iii)  of  Theorem  4.1 

g  At)  =  v(?,0)  e  (I  ),  /?  e  (/?  .+£-1/2, 1/2)  ,  j  =  1,2. 

Q  J.  ^  »  J  l’TJl 

i 

Because  «  >  0  arbitrary  /?  e  (3  .  -1/2, 1/2).  Analogously 

'  J  1+  J “ 1 

for  0  <  *i(/?i  +  1  <  1/2,  g.(i()  e  8?  (I),  '3.  e  ( 3  .  +  .  ^1/2 ,  1 )  . 

*i 

Lemma  4.10.  Let  gJ  e  8*(I),  0  <  3 ^  <  1/2,  0  <  3^  <  1, 


g-  c  8.(1)  ,  1/2  <  13  <  1,  0  <  3-  <  1.  Let  S  =  { r  ,0  1 0  <  Q  <  2w , 

3  1  * 

0  <  r  <  1}  where  (r,0)  are  polar  coordinates  with  respect  to 
(-1,0)  and  $(r)  =  r.  Define 

Ui(r,0)  =  g.(-l+r) 

V.(r,0)  =  «[gi(-l+r)-g.(-l)] 

(by  Lemma  3.1,  3.2,  gi  «  C° ( I ) ,  i  =  1,2,  and  hence  g.(-l)  is 


well  defined) .  Then 


U1  'V1  e  ®5(S>  '  '3  =  '3i  +  1/2 


U2'V2  £  *3{S) '  3  =  V172’ 


Proof .  Assume  first  that  0  <  3^  <  1/2  and  g1  --  8  _  ( 1 ) 
3  =  3^+l/2  and  =  g^(-l+r).  Then  for  k  »  2 


dkU 

.  1,2,  k-2  +  i  . 2  .  _ 

( — £— )  (r  )  rdrdtf 

dr 


<  Cd2k'!g{k)$ 


k-l+*  L2(I) 


_  ,2k.,  ,  ,  2 
s  Cd1  ( k 1 )  . 


Hence  by  Theorem  1.1  of  [4]  we  have  for  k  >  2,  |a|  =  k 


illDa°ll*«+k-21L2(S)  5  °d2k! 


Further 


ii  U  »  <  Ci'cr^1^*  i' 

‘I'lSI  1  1/2  L2(I> 


<  0SaU)i  : 

-  3i  L2(I) 

1 


Hence  U1  e  tt^(S).  Let  now  1/2  <  «1  <  1.  Set  <?  =  1  -  1/2.  As 


before  we  have  for  k  >  2 


2,2.  k-2+/?,2 


( — t—  )  ( r  )  rdrde  <  Cd,  (k!) 


lk(k!  )2 


and  we  get  !:  U  11  <  <® .  Hence  U.  e  35  (S).  Let  us  consider 

*  U A  /  O  >  j- 


H  ( S ) 


now  the  function  V^(r,<?).  Then  as  before 


r  ^  v 

i  /  * 


2,  k-2+^,2 


(  ).rdrd0  s  Cd?k  ( k !  )  2  . 

L  3rk  1 


Further,  using  Lemma  2.5  and  k  >  2  we  get 


4^->2  r-2(rk-2+',)2rdrdO  =  [  l^,2  r'2  ( )  2rdrde 

a a— K“1 


dr  ae 


<  cd2kMg(k-1)*  -1 2 

11  k-2.^  L2(I) 


i  1 1  j  ■  h  j  ■  ■  ■  mjiii  i  in  1 1  |i  m  mi  iji pbjwji  MLPFfjuu  <■■  ■  ipip  p  ■  ip  ii  jornr  [ijr  wii  jyinr  Tnr  rjp  fm 


Lemma  4.11.  Let  g  -=  Std).  0  <  3  <  1/2,  g(-l)  =  0.  then  for 

,3 


0  <  /  <  1/2,  v  =  g$  e  (I). 


Proof .  For  k  i  1 
1 


-1 


((v(k))2*2  -  )dx 

k-l+^+? 


-1 


^=o 

k  .1 


*  _  dx 

k-l+'i+/ 


<  Cd 


2k' 


(g(f)  )2(  (k-€)  !  )2*^_(k_<t)*2  -  dx 


€=0  -1 
k  ,1 


k-l  + i +y 


s  Cd 


2k 


I  <■ 


(f))2*?  (  ( k— f )  !  )  2dx 


*~e  =  l  -1 


?  +l~  1 


2*2  2 
(g)  *.  (k!)  dx 


-1 


4-1 


<  Cd 


2k 


k  ,1 


I 


(g(/))2*-  (  (  k~e  )  !  )  2dx 


^  =  1  J-l 


(g')2*?(k!  )  2dx 


-1 


<  Cd2k(k!)2 


when  we  have  used  Lemma  2.5  in  the  above  inequality.  Further 


r1  r1 

!  2  1  2  2  2 

I  v  dx  =  I  g  ♦  dx  ^  C!'g 

l.j  )_!  Hi'1 


(I) 


3 


by  Lemma  4.1. 


Lemma  4.12. 


Let  g  £  9.  (I).  g(tl)  =  0.  1/2  <  3  <  1,  0  <  /  <  1,2 

3 


v  =  g*  .  Then  for  /}+/  >  1 ,  v  •=•  8 .  (I)  and  for  3+/  <  1, 


-r 

V  £  8?  (  I  )  . 

*+/ 


'3  +/  - 1 


Proof .  (a)  Assume  first  that  3+y  >  1.  Then  for  k  '  2 

.1 


(v(k),2;2 


dx 


-1 


k+/3+/  -2 


<  Cd 


2k 


I 

L-f  =  2 


(g(/> ) 2* 


(  (It-,)  !  )  dx 


-1 


-/-(k~f)+k+3+/  -2 
1 


+  (k!  ) 


g2*?  dx  +  { (k-1 ) ! ) 2 


<  Cd 


2k 


-1 

r  ^  ^ 

k  =  2  -1 


3-2 


-1 


2  ‘  2 

g'^4:  dx 

3-1 


(g(/))2*?  ((k-«)!)dx 

/3+£-2 


+  (k!  ) 


-1 


,  2  r  2  . 

g  ♦ .  dx 
'3-1 


In  the  last  inequality  Lemma  2.5  has  been  used.  Because  by  the 
imbedding  theorem  I g'  ( 0 )  I  *  C  g  -  _  using  Lemma  2.5  once 

H(j  1 1 ) 

more  rendering  that  3-1  >  -1/2  we  get 


!  2*2 

5  g'  ♦ .  dx  s  C 

-!  ^ 


!  g"2*.dx  +  !  g' (0)  !  2 1  4  C1  g  2  < 

3  J  HZ,2(I) 


-1 


fi 


Hence 


(V(k))24 


-1 


k+'3+/  -2 


2k  2 
dx  <  Cdj  (k!  ) 


Further  as  before 


2*2  2*2  2 
v'  ♦.  dx  <  C  g'  4  7  dx  *  C  g  <  <x> . 

*+/-l  J  '?-l  H.(I) 

? 


Because  g  e  C°{i),  v  «e  L2(I). 

(b)  Assume  now  that  3+r  <  1.  Then  for  k  >  2  we  get 


exactly  as  before  that 


Further 


(v(k) ) 2*  .  dx  s  Cd2k(k! ) 2 . 

k+/3+/-2  1 


2  2*2  2*2 
v'  dx  <  C  g  *  ,  dx  +  g'^4^  dx 

-f-1  *  -r 


5  C  g"2*2  .dx  +  Ig'(O)  |2  . 


Because  -f  +  1  >  '3  by  our  assumption  we  see  that 


v'2dx  <  C ii  g !■  2  „ 

,  H.'2(I) 


Using  Lemma  4.2  we  get  also  >i  v  <  C!gi 

L2(I)  H2 ' 2 ( I ) 

Lemma  4.13.  Let  Q  be  a  curvilinear  polygon  with  the  vertices 
Aj ,  i  =  Let  u  e  $2(Q)  and  w  be  such  that 

|Daw|  <  C*_|n|+y |a | !d,a| , 


r  ■  {r 


,/M).  Ia|  e  0,  «1-/1  >  0,  ^  >  0. 


Then  v  =  wu  e  S_(Q)  where  ?  ^  *  ^i~^i' 


a  *  v  *  ■  *  ««  .  *,  •  i  »  _  ^ »  *  ,  *  .  »  _  *  «  ^  ,  *  »  *  •  «»  .  «  ,  «  .  w  »*,  *'  *  „  -  «r_  ./*  «  _  «r_  r,  »■'„  yj  y. 


Proof.  For  k  ?  2,  |a| 


=  k, 


|D|alvlV 


dx  <  Cd 


2k 


Q 


|  a  |  -  2  +>3 


I 


Vr  — ^  ^ 

ID  u|  |D  w| 


^:=o  q 

k 

5  Cd:~  V  ( (*+l)  !  ) 
^  =  0 


*  dx 

k-  2  +  7 


fX 


,_k -f  .2.2  , 

D  u  ♦ .  _  .  dx 

k-2-^  +  < 


0 


Further 


k 

2k  2  2  2k-  2 

S  cdi  2_(('+1)!)  ((k-H-*)!)  <  Cd^  ((k-2 

€  =  0 


| D1v| 2dx  <  C 

I D1u | 2 | w | 2dx  + 

2  12 

1 u |  | D  w |  dx 

Q 

- 

Q 

0 

because  by 

lemma 

2 . 1 

u  e  C°(Q  )  . 

It  is 

very 

easy 

to  prove 

Lemma  4.14. 

Let 

g  e 

*9(1) ,  0  < 

3 

and  v ( i 1 ) 

=  0. 

Let 

g  «  *  -  ( i ) . 

3 

and  v ( : 1 ) 

=  0. 

3  <  1/2.  Then  v  =  g*  8*(I) 

3 

1/2  <  3  <  1  then  v  =  g$  t  8  2 

i 


Proof .  The  statement  that  v  €  **(I)  can  be  directly  verified 

3 

By  Lemma  4.1  v  is  continuous  on  i.  If  v(-l)  *  0  then 
v2(x)  >  e  >  o  for  all  Ix+ll  <  6.  Hence  g2  =  ( v*  *)2  ■  <  $  2  ^ 
which  would  contradict  with  the  assumption  that  g  *?(I), 

0  <  3  <  1/2.  The  proof  of  the  second  part  of  the  lemma  is 
analogous . 


2 

Lemma  4.15.  Let  u  e  8^  (Q )  ,  0  <  tf  <  1  and  u  =  0  at  A  .  Th 
u*  1  **(Q).  The  proof  follows  easily  using  Lemma  2.2. 

Theorem  4.2.  Let  Q  be  a  straight  polygon  with  the  edges  r  , 


$ 

$ 

’.*1 

g 

S 

n[ 

*5 


$ 


i  =  1 . M.  and  let  g  --  *  _  (T  2 )  #  0  <  r  <  1/2,  <  i  =  <.  +  1/2, 

1=1,2  (respectively  g  «?«,),  U%  <<<!.-.  <J-1  2. 


i  =  1,2) 

and  g(Ai) 

=  0,  i  =  1,2. 

Then  there 

is  u  such  that 

(i) 

U  e  *2(Q)  . 

with  0  <  3 

A 

►-* 

C_i. 

II 

3  , 

. . , M ,  arbitrary . 

1 

(  ii) 

uj  =  g 

1 

and  u  j  r  =0 

j 

for  j  = 

2  , 

.  .  ,M. 

i 

M 

Proof.  Let  «  =  |"1 

i  =  3 

l  x-A^  !  2  ,  x  -/  0 

Denote 

g 

=  g ; h  .  Then  obvi - 

1  2 

ously  g  e  $_(r  )  (respectively  g  ® . ( r  )).  Select  now 

3  1  <  1 

0  <  / .  <  1/2  such  that  0  <  <+/ .  <  1/2  (respectively 

1  2  -/  , 


0  <  ■?  +  >  -1  <  1/2).  Denote  g  =  g  ("|  Ix-A. 

1  i  =  l  1 


=  g$  where 


/  *  (^1">2'0'  •  •  •  '0)  By  Lemma  4.1  and  4.2  g(A  )  =  0,  i  =  1,2. 

Using  Lemma  4.11  (and  4.12)  we  see  that  g  -/  (I)  (respec- 

'3+/ 

tlvely  g  e  * l  ( I ) ) . 

3+r- 1 

Let  U  €  H*  (Q  ) ,  AU  =  0  and  U  =  g  on  T  and  U  =  0  on 
f  j ,  j  =  2.....M.  Function  U  exists  and  is  uniquely  determined 
To  see  it  let  <*>(x),  x  e  r  ,  <p  •?  C^r^),  <p  ( x )  =  1  for  x-A^ 
c/2,  i  =  1,2  and  <p(x)  =  0  for  Ix-A^l  >  c,  i  =  1,2  with  >: 
sufficiently  small.  We  define 

U  =  u1+u2 

where  ^  =  0,  UA  e  H1  (Q )  ,  i  =  1,2,  Ux  |f  =  g(l  -<p)  .  U2lr  =  g<r 

and  Ui=0  on  r^,  j  =  2,...,M.  Because  =  g(l-<p)  C  (T 

and  hj(x)  =  0  for  Ix-A^  <  c/2,  Uj  obviously  exists. 

By  Lemma  4.10  there  exists  W  -  H  *  ( 0  )  such  that  W](  =  h2  = 


gp ,  and  W|  =  0,  j  =  2,...,M.  Hence  U.  exists  too.  Function 
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U  has  the  following  properties: 


(  i  )  a  u  =  o  . 

(ii)  uj  =  g ,  u  j  =0,  j  =  2 . M. 

1  j 

(iii)  g  is  analytic  on  r  (not  on  r  ). 

(iv)  in  0.  =  Q  (x||  x-A.  I  <  6  }  ,  i  =  1 , 2  with  o  suffi- 

1.0  1 

ciently  small  there  is  W.  such  that  W.  e  B2(0.  .  ) 

3 

where  3  =  3+r  +  l/2  (respectively  3  =  ?+/ -1  +  1/2)  and 

W. j  -  =  g.  (This  follows  from  Lemma  4.10.) 

1  i  i  ,6 

By  the  selection  of  / ^  we  have  3^  >  1/2,  i  =  1,2.  Using 

now  the  same  arguments  as  in  the  proof  of  Theorem  2.1  in  [4]  we 

2  _ 

conclude  that  U  €  ®_(Q)  where  3.  =  3.  .+1/2  (respectively 

3  111 

'*i  *  '’l^i'1721  '  1  =  1  < 2 '  and  1  >  3^  >  1/2. 

2 

By  Lemma  4.13  we  see  that  u  =  e  3^(0)  where  3^  = 

3  i  +  l/2  (respectively  3  A  =  ^-1/2),  i  =  1,2  and  0  <  i  ^  <  1 

arbitrary  for  j  =  3 . M.  In  addition  uj  =  g  and  uj  =  0, 

1  j 

j  *  2 . M. 

2 

Let  us  outline  the  main  idea  of  the  assertion  that  U  B  (0  ) 


Let  Sj  6  <  r ^  <  6^,  0  <  <  r  ^ )  >'  0  where  (r^.tf^) 

are  the  polar  coordinates  with  the  origin  in  A^.  We  select 
6  <  1  such  that  S.  .  s.  =  0  for  i  *  i.  Using 

0  If  26 ^  j  t  26 j 

Theorems  5.7.1,  5.7.1'  and  6.6.1  of  [17]  we  conclude  similarly, 

M 

2 

as  in  the  proof  of  Theorem  2.1  of  [4],  that  U  -  B  (O  -  u  S.  c 

5  1=1  ’■  1  1 

due  to  the  analyticity  of  g  on  r  -  U  S.  c  , „ . 

1=1  '-V4 

2 

to  prove  only  that  U  Ba(S.  ,  ,  „  )  . 

3  '  i .  /4 


Hence  we  have 


<PQ  e  C  (P  ) 


<PQ  (  1 )  =  1  for  0  i  r  <  1/2 


<pQ(  0)  =  0  for  x  i  1 


<Pb  (r)  =  ‘PqI^ZT'1  =  <p(‘r) 
i  i 


Denote  v  =  <pU .  Then  v  can  be  understood  to  be  defined  on  the 


infinite  sector  =  {{r.,0.)\0  <  r  <  oo,  0  <  0  <<•>.}  when 

W »  11  1  1 


extended  by  zero  outside  of  S.  *  and  we  have  v  e  H  (Q  ) 

1,0i  wi 


Now  we  prove  that  v  ^  R  ( S .  )  as  in  [4]. 

c  /J  l  ,6  . /2  1 


Remark  4.3.  We  have  assumed  that  either  g  e  R_(r,),  o  <  ft  <  1/2 


or  g  e  R_(r  ),  1/2  <  ft  <  1.  Obviously  Theorem  4.2  is  correct  if 

n  * 


1  2 
g  e  *.  (T  )  only  in  the  neighborhood  of  A.  and  g  e  R.(T  )  in 

ft  1  ft  1 


the  neighborhood  of  k^ .  Theorem  4.1  leads  easily  to  the  next 


theorem , 


Theorem  4.3.  Let  Q  be  a  straight  polygon  with  the  edges  P  ^ , 


i  *  1 , . . . , M  and  let 


g  e  *  (P,).  =  (^lrl.«ir2).  0  <  <  1/2. 

i 


Si,l  =  Vl*1'2'  *1.2  ■ 


g  «  *  (Tj) .  S,  -  2).  1/2  <  <  1. 

1 


fll.l  =  8i.l-1/2-  *1,2  *  ,J1.2-1/2’  1  5  «  '=  11 . M>- 


Let  further  g  be  continuous  on  /  =  U  T..  Then  g  t  R_  (>  ) 

ie  Q  1  ft 


i  -  '  i  - 1 , 2  '  i  ,  1 '  '  iur  tti 


arbitrary  for  A^.  e  y  . 


Proof .  Because  g  is  continuous  on  y  we  can  construct  a  poly¬ 


i-1  «r  Q 

or  i  <?  Q 

S 

0 )  and 

0  <  3  .  <  1 

1 

a  poly- 

1 

% 

§ 

construct 

'j 

nomial  P  on  Q  such  that  g-P  =0  at  A^  .  Then  we  can  apply 


Theorem  4.2. 

Remark  4.4.  It  is  obvious  how  the  theorem  may  be  modified  when  g 


1  2 

€  ^  i )  respectively  g  e  in  the  neighborhood  of  A^ 


only.  See  also  Remark  4.3. 


Remark  4.5.  Theorem  4.1  and  Theorem  4.3  are  complementary,  which 


is  analogous  to  the  theorems  of  trace  and  extension  in  usual 


Sobolev  spaces  on  smooth  domain,  namely,  if  g  ?  (T  .  )  ,  0  < 

3  .  1 

l 


3  <  1/2  (respectively  ge  K.  (F  )  ,  1/2  <  3.  .  <  1)  j  =  1,2, 

1  *  J  1 '  J 


then  we  have  an  extension  by  function  G  e  35^(0),  3^  =  3  ^  i  +  i/2, 


'*i  +  1  =  3ii2  +  1/2  (respectively  3  i  =  3^  ^-1/2,  *i  +  1  =  3i  2  - 1  /  2  )  , 


and  if  G  *-i  sj?(Q)  then  G|r  =  g  €  V1.  (T  )  ,  3  -  3  - 1/2,  3 

i  3  .  +c  '  ' 


=  3^  +  ^-i/2  for  1/2  <  <  1  (respectively  g  e  'f> .  ( r  ^ )  , 

3  .  +£ 

1 


'3i,l  =  3  i  ,2  =  fli  +  i  +  1/2  for  0  <  3  i '  3i  +  l  <  i/2)'  f  >  0 


arbitrary . 


Theorem  4.4.  Let  Q  be  a  straight  polygon  with  the  edges  r^, 


i  =  1 . M,  and  let  g  e  *.(r  ),  0  <  3  <  1/2,  i  =  1,2,  ?.  = 

3  1  X  1 


3. +1/2,  i  »  1,2  (respectively  g  €  B.(T  ),  1/2  <  3.  <  1,  3.  = 


3^- 1/2,  i  =  1,2).  Then  there  is  u  such  that 


(i)  ue  $  ^  (Q  )  with  0  <  '?  j  <  1,  j  =  3,...,M  arbitrary. 
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a 


£1 


>; 


>1 


■a 


*1 


I 


I 


I 

*  • 


a 


V  1 
•  .N 

>> 


■Is' 

•  .  I 


(ii)  uj  =  g  and  uj  =0,  j=2,...,M. 

1  j 

"  1  2 

Proof.  By  Lemma  4.14,  g  =  g$  e  ?.(r,)  respectively  _  ( r  ,  )  and 

/3  1  3  1 

--  2 

g(A^)  =0,  i  =  2,3,  and  hence  by  Theorem  4.2  there  is  v  «=:  ft  (0 ) 

such  that  v  =  g  on  r 1  and  v  =  0  on  r^,  j=2,...,M.  By 
Lemma  4.15  the  function  v*  1  has  the  desired  properties.  i 

Theorem  4.4  leads  immediately  to 

Theorem  4.5.  Let  Q  be  a  straight  polygon  with  the  edges  r., 
i  =  1 , . . . , M  and  let 


g  «  *.  (I-,).  0  <  ^i  2  <  1/2, 


*1.1  '  fll,l+1/2'  Si,2  '  ai,2*1/2 


or 


,1 

ii 


g  -  (r,).  ^  ^  i  ,  l '  ^  i  ,  2 ) '  ^  ^  i  ,  1 '  ^  i  ,  2  < 


*1.1  *  "i,r1/2'  ^  i  ,  2  =  ^  i  ,  2 _  1  /  2  '  1  *  Q  ’=  {1 . M> 


1/2 


Let  r  =  U  T  ..  Then  g  e  15  (y  )  where  >i  .  =  max('3.  ,  ■<  .  ) 

ieQ  1  3  1  1-1  ’  ‘  1  '  1 


A^  e  y  (if  i-1  e  Q  or  i  Q  when  we  define  ii  .  _  =0 

respectively  3  ^  1  =  0)  and  0  <  ft  ^  <  1  arbitrary  for  A^  «  / .  a 


Remark  4.6.  It  is  obvious  how  Theorem  4.4  has  to  be  modified  when 


1  2 

9  e  S-(r.)  respectively  g  e  *.(r.)  in  the  neighborhood  of  A. 
3  1  "  1  1 


only.  See  Remark  4.3. 


Theorems  4.3  and  4.5  give  the  characterization  of  the  bound¬ 
ary  conditions  which  guarantees  that  the  solution  of  an  elliptic 
partial  differential  equation  of  second  order  with  analytic  coef- 


ficients  on  a  domain  Q  with  piecewise  analytic  boundary  belong 
to  “B2(Q)  or  p2(Q)  (see  Theorems  3.2  and  3.3). 

In  the  concrete  cases  these  conditions  are  usually  very  easy 

to  check.  Let  us  state  a  useful  lemma  which  characterizes  the 

1  2 

space  $  .  ( I )  (respectively  $  .  ( I )  )  . 

8  8 

Lemma  4.16.  Let 

Qq  =  { z  =  x+iylx  «  I,  lyl  s  a$(x),  a  >  0} 


and  G(z)  be  holomorphic  function  on  such  that  for  i  = 

(t'l't'2) 

I G ( z ) |  <  C*p(Re  z) . 

Let  g(x)  =  Re  G ( z )  J  ^  or  Im  G ( z ) J  j  *  Then  for  >  -l/2+(j-l), 

3  i+P  i  >  1/2+U-D»  0  <  8  L  <  1,  i  »  1,2,  j  =  0,1,2 

g(X)  e  *5(1)  . 

8 


Proof .  By  Cauchy  formula  we  have  for  k  >  0 


(k)  "  "  -k  -k 

g'  '(x)|  <  C*  (x)(*(x))  k !o 


Hence 

1 

*2  „ | g*k* (x) | 2dx 

.  _  k-l+/3 


( Ck  !a 


2 

v  +8  - 1 


dx 


k  2 
(C1dKk! ) 


provided  that  ^  +  8  >  1/2.  Further  for  k  =  0 

|g(x) |  <  C*u (x) 

and  hence  for  >  -1/2,  g  ^  H°(I).  The  lemma  is  proven  for 

j  =  1.  The  proof  of  the  case  j  =  0  is  analogous.  Let  us  con- 

36 


V/ j*  '  j». * 


sider  now  the  case  j  =  2.  We  see  that  for  v  ^  ^  >  3/2  and 


k  >  2 


,  2  .  (k)  .  .  ,  2  .  .  ,  -k,  2 

$  . I g  (x) |  dx  s  (Ck!a  ) 


-1 


k-2  +'? 


-1 


j.--k+k-2+3 


dx  <  (C1dkk! ) 2 


Further  if  v ^  >  1/2  then  also  g  ^  H  (I). 


Instead  of  |G(z)|  i  Cf  (Re  z)  we  can  assume  that 


I G ( z )  -  P(z)|  <  Cf^fRe  z)  where  P(z)  is  a  polynomial 


Lemma  4.16  is  very  useful  in  practice.  For  example  if  g 


is  analytic  on  then  g(x)  can  be  extended  into  some  neigh¬ 


borhood  of  fj,  and  therefore  g  e  5  ?  ( I )  .  Lemma  4.16  characteri 


zes  very  well  the  structure  of  the  spaces  5  (1)  (respectively 


»  3  ( I )  ) 


Lemma  4.17.  Let  g  e  5_(I),  0  <  8^  <  1/2.  Then  there  exists  a 


13 


0  such  that  g  can  be  analytically  extended  onto  Q  and 


a 


G(z)  -  g(-i) (1  x)  -  g(i) (x*1}  1  ^  c* 


.(Re  x) 
1/2 -3 


,0  - 


( g  e  C  ( I )  by  Lemma  3.1). 


Proof .  Since  g  e  5.(1)  we  have  by  Lemma  4.3  for  k 


(i 

,(k) 


(x) |  <  C 


^  -1 


*  -(x) 

^  k-l/2+fl 


V 

d  k! 


Hence  the  series 


0D 


g'(x)  = 


k  1 


^g  ~  *'  (x0)  (x-x0)" -jy,  X. 


,  ( k+1 ) 


e  I 


k=0 


*(x0>  x 


is  absolutely  convergent  for  | x-xQ I 


and  hence  also 


wswu  wawzwawvwvw&e  wvwjwu  pji'w.yjpjr;r.  »■.  rr  rr  r>.  r  rwr  w.  if.1  n'>j  w 


G'  {2 


■I 


( k+l ) ,  w  . k  1 

g  (x0)(z-x0)  ^ 


k=0 


*<V  1  --1 

converges  for  I  z-x  I  <  — ^ —  -  and  |G'(z)l  i  C*.  (x  ),  x  = 

ii +  1/2  U 


Re(z),  find  C  is  independent  of  x^ ,  which  yields  the  lemma. 

So  far  we  have  assumed  that  Q  is  a  straight  polygon.  We 
did  not  exclude  the  case  that  the  internal  angle  is  2 n,  i.e.,  we 
did  not  exclude  the  slip  domain.  Let  us  now  consider  the  curvi¬ 
linear  polygon  and  assume  that  it  is  a  Lipschitzian  domain.  let 
us  prove  first 


Lemma  4.18.  Let  Q  =  (x^x^-l  <  xl  <  1 ,  0  <  x2  <  hfx.^),  h(x1)  > 


a(x1+l),  h(-l)  =  0,  a  >  0).  Assume  that  vlx^Xg)  is  an  analytic 


2  2 

function  on  S  =  ,x2 ! (x^  +  1 )  +  x2  <  4}  such  that 


(i)  rp(xx,  h(x1)  )  =  0. 


dp 


ii)  dx  ( » 0 )  >  ct  >  0 ,  -1  <  x^  <  1. 


Define 


ri  =  {x1>x2!-1  <  xx  <  l,  x2  ■  0} 


r 2  =  (XjjXgl-l  <  x1  <  1 ,  x2  =  h(x1)} 


and  let  T  =  Q  n  where  =  (r,$|0  <  0  <  2 n,  o  <  r  <  1 


where  (r,<?)  are  polar  coordinates  with  respect  to  (-1,0)  and 


*  i 

T  =  S.-T.  Let  g.  e  35  _  ( T  )  ,  0  <  8  .  <  1/2,  13  =  13  (respectively 

1  n  1  1  12 


.  2  ,  r. 


g2  =  35 .  ( r-  x )  ,  1/2  <  I3i  <  1,  81  =  f32),  gi(-l)  =  0,  i  =  1,2  and 


f  =  r . 

Then  there  exists 
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■  ’  J>  *>  *•  ’>  *_*»  *Ji  '  >  r  J>  *>  V#  *.-**•  *  J*  *j«  *  m  *  m  .  *  .1 


,**  ,  *  d  » 


•J’  X->  'JTJV  fW  *5*  'WW  rWlRTW. 


'wiwwwwn »  upu  pi* 


V  e  B2(T),  V*  e  B2(T*),  3  =  3  +  1/2 

3  3 


(respectively  V  e  B  (T),  V.  e  B  (T  ),  3  =  3-1/2) 


3 


3 


such  that  Vi  =  and  =  gA  on  r  o  t  and  Vi'Vj_  =  0  on 


r  2  °  T. 


Proof .  Let  y>(r,<9)  =  v{r,0)  — .  Then  <p  ( r,0)  =  ^(x^)  is  analytic 


-1 


on  and  tp(K^)  >  a  >  0,  hence  <p  (x^)  is  analytic  on  f1 


a 


too.  In  addition  <p  =  0  on  f  2>  Further  |D  <p(x.,x2) 


C  jo  !  !  $  l°lcjlal  by  Cauchy  theorem  of  the  theory  of  two  complex 


^  - 1  ~  1 

variables.  Define  gx  =  g ^  (x1).  Then  gx  €  « r (f  x )  and  by 


Lemma  4.11  there  exists  U  on  S  such  that  U  e  B  (S 

•  •  ^  -T 


'3  +  1/2  and  U^Jp  =  •  Define  now  V.^  =  U.^ .  Using  Lemma  4.13 


we  conclude  that  V  e  B2(T)  (respectively  B2(T*)),  V  |  =  g 

3  3  1  1  1 


and  V 


•r  =0.  The  proof  that  V  has  desired  properties  is 

^  1  O  & 


quite  analogous. 


Lemma  4.19.  Let  Q  =  (x^x^-l  <  x:  <  1 ,  h1(x1)  <  x2  <  h2(x1 


h^Xj)  <  -atXj  +  l),  h2(x1)  >  a(x1  +  l),  a  >  0,  h.(-l)  =  0,  hi(x 
analytic  functions  on  1,  i  =  1,2}  and 


ri  =  (x1,x2!-1  <  x1  <  1 ,  x2  =  hi(x1)} 


=  D  o  s  ,  =  (r,0|O  <  6  <  2-n  ,  0  <  r  <  3  ,  3  >  0}, 


0  =  S  -0  , 

3  3  3 


where  (r,3)  are  polar  coordinates  with  the  origin  at  (-1,0). 
rl 


Let  g1  e  B*(T  ),  0  <  3  <  1/2  (respectively  g  ^  B2([  ), 

3  1  2  3  1 
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Li ru 'rvnrv»TJiruim,"ir\i w\) /vri  w -v  wn/rw 


1/2  <  «i 

<  1  )  ,  S 1  =  8  2  ,  g 

exists  r? 

>  0  and  e  ! 

arbitrary 

( respectively 

such  that 

VilrioQ??  =  9i 

Proof.  Because  h  (Xj) 

extended 

onto  I§  =  (-1-6 

(x1,x2)  — 

(Yi'Ya*'  yi  "  X1 

t?  =  6/2 

and  M(Q  )  ■  Q  . 

T)  T) 

*  * 

<->  s  =  r 

> j  T  ^  where  f 

** 

Then  there 


fl  77 


fl  77 


*  77 


B  77 


311(1  Vi  If  nO  -  °‘ 
2  7? 


sufficiently  small  we  have 


dQ 


1  =  iY1'Y2\-l  <  Y1  <  -1+7?,  Y 2  =  °>< 

J. 

r  2  =  {yj-Yg!-1  <  Yi  <  -1+^!'  Y2  =  h*^)  =  h2(Y1)-h1(Y1)  and 

h2^Yl^  >  °,1(Y1  +  1)*  In  addition  it  is  easy  to  see  that  |J|,|J  1| 

<  //  <  oc  where  J  is  the  Jacobian  of  the  mapping  M.  Because 

*  _ 
h2(Yi)  is  analytic  on  -l  <  y1  <  -l+rj^  we  define  «’(y1,y2)  = 

* 

-y2  +  h2^yl^  and  'My1(y2)  has  the  properties  in  Lemma  4.18. 

Using  now  Corollary  4.4,  4.5,  g  e  <r.(T  ),  g  e  (r  ?  ( r  )  and  hence 

8  1  8  1 

using  Lemma  4.6,  g  ( M~ 1 ( y ) ) j  e  <r*(T*),  g  (M_1(y 


■)  > 


y2=o 


2  * 

®-(r.).  Using  Lemma  4.8  and  Lemma  4.18  there  are  functions  V. 

a  1 

*  * 

and  (respectively  V2  and  )  on  n  (respectively 


S  )  (respectively 
7?  o 


Q  °  S  ) ,  which  belongs  to  R.  (0 

7  V  2  8+c  /  2  77  c 

2  .  * 

®-  (v  n  S  )).  Using  now  Lemmas  4.7,  2.3,  our  lemma  follows.: 
(l+c/2  1  71 2 

The  lemma  leads  to  the  following. 


Theorem  4.6.  Theorems  4.3  and  4.5  hold  also  for  Lipschitzian  cur¬ 
vilinear  polygon  when  8^  are  replaced  by  /Jq+jt  ,  f  >  0  arbitrary 

Proof .  Because  the  edges  are  analytic  curves  and  g  are  analytic 


-.v.v.v.v.v.v.vv 


\-WV. 
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s 


on  r0  (but  not  on  fQ)  we  show  similarly  (as  in  the  proof  of 
Theorem  4.1)  that  the  solution  u  of  the  Laplace  equation  belongs 

to  $2  (Q  )  .  This  can  be  done  identically  as  in  the  proofs  of 

3+t 

2 

S+£ 


Theorems  3.3  and  3.4  of  [6],  showing  that  ff “  _ (Q ) • 


Remark  4.7.  Comparing  the  respective  theorems  for  straight  and 
curvilinear  polygons  we  see  that  in  the  latter  case  we  are  losing 
slight  in  the  regularity.  It  is  not  known  whether  this  loss  can 
be  removed. 


5 .  The  finite  element  method 

Let  us  consider  the  finite  element  method  for  solving  the 
model  problem  (3.1).  We  will  assume  that  Q  is  a  curvilinear 
polygon  and  for  simplicity  of  the  exposition  we  shall  assume  that 
the  vertex  A.^  is  located  in  the  origin  and  the  singularity  occurs 
only  in  the  neighborhood  of  A^.  In  this  case  we  can  assume  that 
#  =  r . 

Let  us  first  describe  the  meshes  which  we  will  consider.  Let 

Qn  =  {« •  ,,  j  =  1 , . . . , n+1 ,  i  =  l,....I(j))  be  the  partition  of  0 

1  >  3 

satisfying  the  following  conditions  (see  Figure  5.1  where  indices 


i,j  of  Q.  .  are  given): 

1  t  3 

(i)  Q.  are  open  quadrilaterals  or  triangles  (curvilinear 
i »  J 

quadrilaterals  or  triangles),  the  intersection  of  any  two  CK  j 
is  a  common  vertex  or  the  entire  side  or  is  empty  (the  mesh  shown 


in  Figure  5.1  is  a  geometric  mesh  with  respect  to  the  vertex  A^^ 
(see  (iv)).  If  the  singularity  would  occur  also  in  other  vertices 


then  similar  refinement  would  be  in  the  neighborhood  of  A ^ ,  j  > 

1)  . 

(ii)  Let  h,  .  be  the  diameter  of  Q .  .  and  h.  .  the 

*  t  J  *  '  J  1  •  J 

diameter  of  the  largest  circle  inscribed  in  !) .  .  .  We  shall 

*  '  J 

assume  that  there  is  a  constant  \  independent  of  n  such  that 
(5.11 

(ill)  Let  M  =  (M  ,  1  5  1  *  1(3).  i  s  3  s  ii  which 

^  *  3 

Mj  j  is  one  to  one  mapping  of  the  standard  (master)  square  S  = 

[1,1]  '[-1,1]  respectively  standard  triangle  T  =  {K,ri\  O'  ri  ' 

1-? ,  -1  <  f  <  1)  onto  0.  ..  If  T  is  a  triangle  then  we  will 

i  •  3 

assume  that  M.  .  can  be  extended  into  standard  square  S  (T  is 

*  •  J 
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Figure  5.1.  Scheme  of  the  mesh. 


half  of  S)  such  that  M.  .  (S)  «  G.  .  c  Q  and  M.  .  still  sat- 

*  <  J  ^ « j  1 1 J 

isfies  on  G^  j  all  conditions  which  will  be  later  imposed  on 

Mi  j’  Let  pj  j  (  and  -y  i  j  y  denote  the  vertices  and  sides  of 

Qi,r  then  Mi!j(pi,  and  Mi!j{^i,j,^  are  vertices  and 

sides  of  S,  1  ^  l  <  4  (respectively  vertices  and  sides  of  T 

with  1  <  f.  s  3).  Moreover  if  j  and  Mm  ^  maP  (closed) 

standard  square  S  onto  element  Q,  ,  and  Q  .  with  the  common 

l  #  j  m ,  k 

side  <t>  -  P  ,P  then  for  any  P  e  y ,  dist  (M_1  (P)  .m”1  (P  )  )  = 

^  *•  i  <  j  i  i  J  * 

t  =  1.2. 

Let  M^  j ( S )  Mffl  ^(T)  =  r  ~  be  a  common  side  of  the 

quadrilateral  and  a  triangle.  If  y  is  the  image  of  the  sides 


where  c.c.v  i  =  are  positive  constants  independent  of 

i,j,n.  If  O  is  a  triangle  then  we  assume  that 

*  •  J 

(5.9)  dist { G .  .,0)  -  C  dist(0.  .,0) 

i  <  J  i  -  J 


and  if  /  is  the  common  side  of  Q  .  and  0  ,  j  >  >,  then 

i  ,  j  m  ,  * 

is  the  side  of  G.  .;  if  y  is  the  common  side  of  <>  and 

1  <  J  '  .  ] 

Q.  .  ,  then  y  is  the  side  of  G.  or  G  if  /  is  the 

1  >  J  l  .  J  *  >  J 

side  of  Q.  .  and  the  part  of  <K>  ,  then  }  is  the  side  of 
1  •  J 

G.  ,.  In  Figure  5.2  we  show  the  association  of  G.  and  . 
i. j  i .  J  i  ■  J 

In  our  example  R ( 1 )  =  5 ,  R ( 2 )  =  1 2  as  can  be  seen  the  numbering 


is  largely  arbitrary.  ^  are  shadowed  by  full  lines  and  G ^ 

by  extended  dashed  lines.  The  indices  i,j  are  indicated  in 
Figure  5.2. 


Let  us  verify  now  our  assumptions.  The  condition  (5.9)  is 


obviously  satisfied.  Let  y.  -  Q,,  .  0.  Then  is  the 

1  4,1  1 

side  of  G^  2  (=  2  'J  which  is  our  condition. 

Let  =  j  °  ^2  l’  Then  f  2  neither  side  of  1 

nor  G2  j  and  our  assumption  is  not  satisfied.  In  this  case  we 
have  to  define  G.  =  Q  u  0 

*/l  *  /  1  4 2 

In  application  we  can  always  assume  that  a  proper  association 
between  ^  and  ^  always  exists.  Nevertheless  we  remark 

that  our  assumptions  mentioned  above  could  be  difficult  to  precise 


ly  verify,  especially  that 


is  one  to  one  mapping. 


theless  this  is  common  in  the  finite  element  practice. 


Never - 


Figure  5.2.  The  scheme  of  the  mesh,  Q.  .  and  G.  .. 

1-3  i<3 

So  far  we  have  assumed  that  the  geometric  mesh  was  refined 
only  in  the  neighborhood  of  one  vertex  (singular  point).  Analo¬ 
gously  we  define  the  geometric  mesh  in  the  neighborhood  of  every 
or  some  vertices.  Instead  of  the  formal  definition  we  show  in 
Figure  5.4  a  geometric  mesh  for  the  domain  Q  shown  in  Figure  5.3 
The  vertices  in  which  neighborhood  the  mesh  should  be  refined  are 
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numbered . 


Figure  5.3.  The  domain  Q. 


Figure  5.4.  Geometric  mesh  on  the  domain  <>  shown  in  Figure  5.3. 
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Let  now 


P={p.  .  ,  1  <  i  i  R(j),  1  '  j  n+1  >  and  Q  = 

^  •  J 


;q.  .  ,  1  <  i  <  R ( j ) ,  1  «  j  <  n+1)  be  the  degree  vectors  with 

1  <  3 


integers  ^.q^  j  ?  We  the  subspace  S^’^(O^)  = 


.-1 


(<p\<p(x1,x2)  =  §(Mi  j(x1<x2))  for  (x1#x2)  €  0.  ^  ^  ,  /? )  ,  K  ,  n 


e  S  (respectively  Z  ,1  e  T)  is  a  polynomial  of  degree  <  p. 

1 1  1 


in  F  and  of  degree  ?  q.  .  in  r?  ,  (respectively  of  total 

t » 1 


degree  max(p.  .  ,q.  .)).  Further  we  denote  s— ' ^ ' 1 ( o n )  =  S-'^(0n; 

•l  #  j  i  #  3  o  o 


°  H  (Q)  (usually  but  not  always  p.  .  =  q.  . 

1  f  J  1  •  J 


Let  us  impose  now  additional  assumptions  on  0 


n 


a 


First  let 


n 


us  assume  that  0 .  .  e  Q 

i ,  3  ° 


are  quadrilaterals .  In  this  case  let 


,  n 


Y  ■  i  o,  l*-fs4  be  the  side  of  the  quadrilateral  0  .  .  e  0  . 

l ,  3  ,  <■  l  ,  j  o 


Then  we  assume 


(5.9a) 


1.3.* 


x  =  h  <p  .(Z) 

**  /  J  ^  *  J  f  ^ 


Y  =  hi.3*’i.3.<(,:1' 


-1  s  t  s  1  ,  (  »  1 , 3 


(  5 . 9b ) 


-1  s  7?  s  1,  (  =  2,4 


and  that  for  some  constants  C  ?  1,  L  >  1,  which  are  independent 
of  i,j,^  we  have 


(5.10) 


1 9 


(k; 


(k) 


<  CL  k! ,  k  =  1,2,... 


i.j.*"  i,3>*' 

and  that  the  mapping  M.  .  which  maps  S  onto  Q  ■  H  has  the 

i »  j  1 '  J 

form 
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Figure  5.5.  The  curvilinear  quadrilateral  0 
and  the  standard  square  S. 
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In  the  case  that  0.  .  is  a  triangle  the  mapping  is  essentially 

*  >  J 

similar.  We  will  define  it  only  for  the  case  when  only  one  side 
is  curvilinear.  Figure  5.6.  shows  the  notation. 


i.j.3 


ili.i 


l.i.2 


Pi.j.2 


r,  p2 


Figure  5.6.  The  curvilinear  triangle  Q 

1  t  J 

and  the  standard  triangle  T. 


X  =  X  (€  ,T)  ) 
1  '  J 


(5.11)  M.  .  = 
1  <  J 


y  -  v 


f  f  .  1-t  1+0  1 

=  ((♦(fl-Xj  —  -  x2  —  .} 

+  ♦  x2  lii  +  x3  § 

=  {(*<'= '-»!  ¥  -  T^Wh.J 

fl-^-O  K+l  J7 

+  Yl[— 2—  J  +  Y2  -r*  +  Y3  7’ 


We  see  that  we  can  extend  M.  .  onto  the  standard  square  S. 

*  '  J 

Let  us  now  describe  the  finite  element  method.  It  is  a  stan 
dard  one. 

(a)  First  given  the  nonhomogeneous  Dirichlet  (essential) 
boundary  condition  g^°^  on  we  project  it  into  the  space 

on  traces  of  the  subspace  SP ' ^ ' 1 ( 0^ )  =  S.  We  denote  this 


projection  by  gg ,  i.e.  we  replace  g 

(b)  The  finite  element  solution  ue 


on  f  by  gc 


S-'Q' 1 (0nJ 


is  now 


defined  in  the  usual  way  such  that 


B(u_,v)  =  fvdx  + 

o 


[1]  , 
g  vdx 


holds  for  all  v  e  S-'^,:1(Og)  H^(Q)  when  ug  =  gS  on 


(5.12)  B(us.v) 


f  V  s  dv  A  Vk  us 

[  2-  3i  .  j  dx .  dx.  2-.  ^i  dx.  V  + 
i . j=l  J  i=l 


cuv  dx 


We  are  assuming  that  B(u,v)  satisfies  the  usual  inf  sup  ( B — B ; 
condition  (with  the  positive  constant  independent  of  SP'^'1(Q^; 
on  hJ(Q)-hJ(Q). 


The  projection  g 


gg  is  possible  to  define  in  different 


ways.  See  [10],  [11].  We  will  use  the  projection  analyzed  in  [9 


Let  y  =  M .  Ay)  c  r 
>  J 


where  y  =  (-1  <  K  <  1,  rj  =  0)  and  let 


g{f)  -  g(M  .(/)).  Then  g  is  defined  on  /  .  (Because  we  will 
*  *  J 

assume  that  g^°^  e  $*(/),  ft  <  1/2,  it  is  continuous.)  We  now 
define 

Sfs  ( ^  )  =  a+bf  +qp 

where  qp  is  a  polynomial  of  degree  p  =  p^  ^  in  K,  qp(-l)  = 
q  (1)  =  0  and  g_(£)  =  g(£)  for  £  =  ±1.  Polynomial  q  is  now 


such  that 


(5.13 


P-1 

«p  =  IVk'^1 

k=l 


where  ^(J!)  are  Legendre  polynomials  and  b^  are  the  coeffi- 


cients  of  the  Legendre  expansion  of  (g-a-b<)'.  We  mention  that 
the  sum  in  (5.12)  starts  with  k  =  1  because 


1 

(  g-a-bK  )  'dx  =  0  . 

-1 

Further  we  underline  that  for  any  g{K)  which  is  continuous  we 

define  b^  in  (5.13)  by  the  integration  by  parts. 

Finally  we  define  gc(x)  =  g* { M . 1  .  ( x ) )  .  We  have  now  the  fol- 

b  b  1  ,  J 

lowing . 

Theorem  5.1.  Let  Q  be  a  polygon  or  curved  polygon.  Assume  that 
the  solution  u  of  problem  3.1  belongs  to  <S^(Q),  ,  g^°^  is  con¬ 

tinuous  on  r  t0]  and  g.  =g^°^/T.  e  35  _  ( T  .  )  ,  0  <  3  .  <  1/2  or  g. 

l  i  %  i  1  1 

e  35  ?  ( T  )  ,  1/2  <  3.  <  1,  i  e  Q.  Let  S  =  SP,Q,1(Qn)  and  u  is 
(i  1  1  ab 

the  finite  element  solution  defined  above.  Let  0  <  u  <  v  <  & 

and  let  *0  s  PA  ^  =  qi ,  j  <  l>n'  pi,j'qi,j  i  1  •  Then 

-hN1/3 

(5.13)  I!  u— u_  I'  .  <  Ce 

H1  (0  ) 

where  N  =  dim  S  '^'*(0^)  and  the  constant  C  in  (5.13)  is  inde¬ 
pendent  of  N. 

The  proof  will  be  given  in  the  next  section. 

In  Theorem  5.1  we  assumed  that  the  solution  is  in  the  space 

2 

&q(Q).  In  section  4  we  discussed  the  structure  of  the  input  data 

in  the  boundary  condition  functions  g^"^,  i  =  0,1.  g^^  and  f 

guarantee  that  the  solution  of  problem  (3.1)  belongs  to  the  space 

«>• 

Remark  5.1.  Assume  that  b^  =  0  and  c  >  0  in  (5.12).  If 


gt0-*  =  0  and  SP2,Q2,1{On)  r>  SPl'Ql'1(On)  then 

a  a 

(5.14)  B(u  ,u  )  <  B(u  , u  ) 

bl  bl  b2  b2 

where  u<,  e  SPl  '  ,:1(Qn)  Is  the  finite  element  solution.  If 

1  ' 

g^°^  *  0  then  (5.14)  does  not  hold  in  general.  If  g^1^  =  0  and 

f  =  0  and  g^  belongs  to  the  space  of  traces  of  SPl  '  '*(0^)  , 

i  =  1,2,  (i.e,  g  =  g  )  then 

S1  S2 

(5.15)  B(u  ,u  )  >  B(u  ,u  ). 

bl  bl  b  2  b2 

If  g^0^  does  not  belong  to  the  space  of  both  SPl  '  (^i  '*(0n),  i  = 
1,2,  then  (5.15)  does  not  hold  in  general  (although  numerical 
experience  shows  that  in  most  cases  (5.15)  still  holds). 


6. 


The  rate  of  convergence 


We  will  prove  in  this  section  the  statement  of  Theorem  5.1 
Let  us  prove  some  auxiliary  lemmas. 

Lemma  6.1.  Let  g  e  Hk(I),  k  >  1,  g(-l)  =  g ( 1 )  =  0  and  let  p 
2 .  Let 


(6.1) 


g(x)  =  Xa/: 


where  -r(x),  j 


=  0,1,...  is  the  Legendre  polynomial.  Let 


(6.2) 


^ _ * 

f '  ( x )  =  a  .£  . 

p  Z.  j  j 


(6.3) 


g  (x)  =  g ' ( x ) dx . 
P  P 


Then 


(6.4) 


gp(±  1 )  "  0 


(6.5) 


lg(m,_g(m)  ||  2 

9  9p  L2(I) 


<  r  1  r (p-s+l), _(s+l)  2 
-  p2  (  1-m)  f  (p+s  +  1  )  '•  9  L 2  ( I 


p  >  Si  0,  m=0,l 


Proof .  Because  g'  e  L2(I)  expansion  (6.2)  exists  and  because 
g(±l)  =0,  aQ  =  0  and  (6.4)  holds.  Further  obviously 


lg'-3p"L,  *  X'V 


j=p 


We  have 


/..x  i  r  ( j-s+i)  ns,s,.., 


where  (x)  is  a  Jacobi  polynomial  with 


(l-x2)sP®'s(x)P®'s(x)dx  = 


0  for  j  *  k 

s2s+1r  2 ( s+ j+i ) _ 

(2s+2j+i)r ( j+i)r ( 2s+ j +1 ) 

for  j  =  k 


Hence  we  get 


which  yields 


(s+1)  ,2 


2  =  .i  \  a  /  ^  s  ^  (  x )  ;l  2 

l2(I)  Z-  j  j  (  ,il2{I) 


fl  r  ®  ■>  2 

. ,  2.  s  V"  „(s)  .  .  , 

>  ( i-x  )  2Lai3  x  dx 

-1  *-j=s 

>  V*  |  a  I  2  2  ^  ^  • 

'  Z-  1  j‘  2  j  +  1  (j-s)! 

j=P 


la'-a'll2  <  (p~s)  !  Via  l  2  2  i  ! 
ig  9P  L  ( I )  -  (p+s)!  Z-  j  2 j  +  1  (j-s)! 

j=P 


<  l£l£li||gS+1l|2 
-  ( p+s )  !  9  1  L_  (  I 


and  we  get  (6.5)  for  m  =  1.  Further  we  have 


3_3p  = 


which  immediately  leads  to  (6.5)  with  m  =  0. 

Lemma  6.2.  Let  g  e  H  *  '  1  { I )  ,  0  <  3  <  1/2  (respectively  g 
H2 ' 2  ( I )  ,  1/2  </?<!)  an  g  be  defined  by  (6.2).  Then 


C-fgp'!  g;! 


JP  l2(D 


h*,;L(I) 


respectively 


,9P’1,(I)  1  CpJg"„2,2 


h:  •  *  < 1  > 


!Vv»  '  “9p!9!!H^m 


Proof .  We  have  for  o  <  8  <1/2  and  g  e  ' 1 ( I ) 


a. |  *  | 

j  2 


g'  (x)f  j  (x)dx| 


-1 


<  111!  ng:i 

2  ' 1  ( I ) 

<  C^|i||gll  j 

'  1 ( I ) 


Kj (x) | 2*~2 (x)dx 


-1 


because  Kj(x)|  <  1.  Let  now  g  e  H2 ' 2 ( I ) ,  1/2  <  8  <  i. 
also 


.  2  j  +  1 

3  j  5  "f- 


g'  (x)-f  j  (x)dx| 


-1 


*  c llil  .i 


2  l9"„l 


„2  j  +  1,. 

«  c — 5 — 11  g1-  7  , 

H‘(I)  ^  Hg*  (! 


Hence 


,  >  ii  2 


P-1 


JP  1*2  ( I ) 


X  lak'2  TE+T  *  Cp2,,lgii^l.i 


k=l 


H«  (I) 


(respectively  ( Cp2 >! gl! 2  5  ),  and 

H^'d) 


lBp'v*>  5  cI'VaT^’  0Kw,*i,b,b2-*(I 


1/2 


Then 


Proof  of  Theorem  5.1.  The  basic  idea  of  the  proof  is  very  similar 
to  the  proof  of  Theorem  5.3  of  [6].  Hence  we  will  outline  only 
the  basic  steps  and  underline  the  essential  differences  in  detail. 
For  simplicity  and  without  any  loss  of  generality  we  shall  assume 
that  there  is  a  singularity  in  one  vertex  A1  only  which  is 
placed  in  the  origin;  we  did  make  the  same  assumption  also  in  [6]. 

We  shall  first  assume  that  the  mesh  consists  only  of  the  quad 
rilaterals  and  that  p.  .=q.  .  =  p.  >  1.  The  proof  has  few 

steps  similar  to  those  in  [6]  ,[13]. 


Step  1.  Denote  U.  ..  (?,r?)  =  u(M.  Then 

1  •  J  1 1  3 

(i)  for  j  >  1,  U,  .  is  analytic  on  S. 

i »  J 

(ii)  for  some  constants  d  and  c  independent  of  i,j,  i  = 

j  =  l,2,...,n+l,  and  |a|=k,  k=l,2,...  we  have 


(6.6) 


The  proof  is  given  in  Lemma  5.1  of  [6]. 

Let  r t  -  1 ,  .  .  .  ,  4  be  the  sides  of  S.  Assume  that  ^ 
lies  on  K  axis  (i.e.,  =  I)  and  v(K)  is  a  polynomial  of 

degree  p  on  y ^  and  vanishes  at  the  end  points  of  y^.  Then 
there  exists  polynomial  V(^,r?)  of  degree  p  in  K  and  r/  such 
that 

(6.7a)  V  ( *(  ,  17  )  |  =  v(K  )  ,  V  (K  ,tj  )  |  =  0,  l  =  2,3,4 

'  1  r  l 

and 


(6.7b) 


'I  VII  .  <  C|v|  , 

H 1 ( S )  n  (y 


For  proof  see  Lemma  5.2  of  [6]  or  [13]. 


ritmKKXTTTDnCTV FTHm'JAUJU*  J' 


Step  2.  We  construct  polynomial  4  .  ..(?,)?)  of  degree  p,.  in  *■ 

i ,  j  j 


and  n  on  S  such  that  (.  .  =  #.  .  at  the  vertices  of  S  and 

1  >  J  i ,  J 


for  m  =  0,1,2,  1  s  s^  s  p  . 


_  „  9  (p.-s  )! 

(6.8a)  |!D  (U-  .-*  .  . )  I1  _  5  C  - - * - —  i;  U  .  .!: 

t.i  1/J  H°(S)  P j+S j  +  2“2m)  !  1.J  HSJ+3 


2 

"'(S) 


and  using  (6.6) 


<  c 


f  .  r  [dSJ  t3(s  .«>  !]  V<  1~’  »  (n-j  +  2) 

(p  +s.+2-2m)!  1  j  1 


J  J 


(for  the  proof  see  Lemma  4.2  of  [13]).  Define  <o  .  ,(x,y)  = 

i  »  J 


.-1 


for  J  !  2: 


then  we  have  for  0  <  m  <  1 


H  Dm  <  III  <  Ch]'”:iDm(U  -*  ):l 

1,3  H  (Q.  .)  ' 3  ' 3  1,3  H  (s) 

*  •  J 


?hi-»r  (p.rsj|!  i,/2dB,+3  +31,ctu 

Chi. jlcPj+Sj+i-Zm) ! J  d  1  3  31 


~B) (n+j  +  2) 


For  j  =  1  we  use  p  =  1  and  get 


(6.9a) 


,  „  1  1-8.. 

i,1  Uf  #  ,  ^  5  Ch^  i ,:U':  2  2 

,  '  H  (S)  1,1  H^'  (0 ) 


(6.9b) 


^  i  j.  *  j_  i  ^  '•  2  2 

1,1  1,1  H^(S) 


Ci!  ul! 


2  2 

H/J  '  Q ) 


(6.9c) 


))  u-?o  j  .  J!  , 

A  ,  A 


s  Ch}  1  )!u!;  ,  7 
H"(0i  x)  1,1  H‘*^(0) 


,2,2 


in  (6.9b)  we  define  '  (S)  with  the  weight  $  =  r  with  respect 


to  (-1,-1)  and  we  assume  that  M.  ((-1,1))  =  (0,1)  =  A 

■1  /  «!■ 


1  ‘ 


Step  3.  The  function  <p  t  are  constructed  separately  on  every 

^  <  3 


Q,  .  (hence  the  function  <p  composed  from  <p ,  .  «  H  (Q )  )  .  Let 

i ,  j  1 ,  J 
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us  assume  now  that  /  =  Q.  .  ■  Q.  Then  (<p .  .-</>.  .)  |  =  »;  *  0 ; 

J-  #  J  ^  i  »  J  ^  '  / 

nevertheless  v  =0  In  the  end  points  of  }  .  Let  us  assume  first 
that  j  >  €  >  2  and  that  y  =  NL  ^(1).  Denote  (»('•)  =  'i  { NT 
Then  y>[K)  is  a  polynomial  of  degree  p  =  max(pj,p^,)  on  I 
vanishing  at  ± 1 .  Using  the  imbedding  theorem  we  get 


"  V  ••  1  £  C  max(-:U  -*  .  ,  U  -4  ) 

H  (I)  1,J  H  (S)  '  ’  H  (S) 

and  hence  by  (6.7)  there  is  a  polynomial  V(^,r?)  of  degree  p 

such  that  V.  .  (?,r?)  =  i>  on  I  and 
i ,  j 


i!V  U,n)\: 

H  (S) 


H1  (  I 


We  estimate  then  by  ( 6 . 8 )  for  m  =  2.  For  j  =  *'  =  1 

H1(I) 

function  »/>  =  0  on  y  .  For  j  =  2  and  *-  -  1  we  proceed  simi¬ 
larly  using  (6.9b).  In  this  way  we  construct  correction  function 
VA  and  VA  j(x1,x2)  =  V.  j  (  mT  1  j  ( x  x  ,  Xg  )  )  so  that 


<P .  ,+V .  are  continuous  on  every  r  *  T,  r  -  0.  . 

i .  J  -l »  J  ^  <  J 


i  .  e  .  , 


the  composed  function  <p  such  that  <p 


_  =  <p  .  .  +V .  .  belong  to 

,  1.3  1.3 

»  J 


Ha(Q)  and 


R(l) 


I  u~v> !:  2 !  s  c[l!u!;2  ]T  h2(J' 

HX(0)  L  Ht'NO)  f—  1'A 


n+1  R ( 1 ) 

I  I 

j=2  i=l 


,lprVI,  (dsJ+3(S.*3).i2(.2(i-,,,nt2- 

(  P j+s j~2 )  !  J 


Step  4.  We  estimate  now  u-<p  at  the  boundary  r  Let  /  = 


Q.  .or'  .  Assume  first  that  j  *  2.  Then  using  (6.5)  and  the 
I »  J 


assumption  about  g 


we  can  construct  in  the  similar  way  as 


59 


so  that  (affix  the  correc- 


before  the  correction  function  V.  . 

*  •  J 


tion)  =  g  on  r  ^ 0  ^  -f  ^ 0  ^  U  0.  and  (6.10)  still  holds 

j=l  1/J 


1 f  '  "  °i.l  '  r 


[0] 


then  we  use  Lemma  6.2  and  the  assumption 


that  p  s  in  and  analogously  as  before  we  construct  the  correc¬ 


tion  function  V.  so  that  function  <p  H  (0  )  is  constructed 

l  /  l 


which  has  the  following  properties: 


(i) 

<P 

e  HP ' 

Q,  l 

<  ii ) 

<P 

=  gs 

on 

(  iii! 

l 

ii  u-<*>l 

.  2 
'.,1 

< 

/  /A  \ 

C  n 

<°;>- 


(0) 


2  2  (  l-'?  )n 


(6.11) 


n+1  ,  .  . 


j~2 


(Pj+Sj-2) 


iSj+3,_  +  g  ,  j  ,2tT2(1-J)  (n+2-j) 


-  ( d  J  "(S.  +  3)  !  ) 


In  (6.11)  we  have  used  the  assumption  about  the  mesh,  namely  that 
R(j)  <  K  independently  of  n. 


Step  5.  So  far  we  have  not  chosen  in  (6.11)  the  values  of  s^ . 


By  the  same  procedure  as  in  [6]  we  can  select  s ^  in  dependence 


on  Pj  so  that 


(6.12) 


\i-f  <  C  e 

H1 (0  )  1 


-bn 


2  3 

and  because  N  <  K(nn)  n  <  Kjn  we  get 


(6.13) 


~y  N 

?  C  e 

H  (0) 


1/3 


Step  6.  Let  now  u  be  the  exact  solution  of  problem  (3.1)  such 

,  C  o  ] 


that  g' 


is  replaced  by  gs*  Then  ug  =  u-u  satisfies 


60 


•  %  *.  V  N  %« 


As  in  Step  4  we  construct  function  v  e  H  (0 )  such  that 


and 


[0]  _(0) 
v  =  gs-gL  on  r 


V"  <  C"u-<p 

H1(Q)  H  (0 


Because  we  have  assumed  that  the  bilinear  form  associated  to  prob¬ 
lem  (3.1)  satisfies  the  inf-sup  (B-B)  condition  we  get  immedi¬ 
ately 


(6.14) 


(,'u  >!  s  C))  u-<o!l  . 

S  HA(Q)  ir(Q) 


Step  7.  Finite  element  solution  u0  can  now  be  understood  as  the 
finite  element  solution  of  the  problem  with  exact  solution 
u  =  u  +  u_ .  Now  we  have  using  (6.13)  and  (6.14) 

w 

«•  1  /  3 

H  u  If  <  Ce  ' 

5  H1(Q) 


and  hence  also 


Therefore 


l!u-u  "  .  <  Ce 

55  H1  (Q  ) 


~7  N 


1/3 


u-u  .  <  Ce 

H  (ft  ) 


-y  N 


1/3 


which  was  to  prove. 


So  far  we  have  assumed  that  the  mesh  consists  only  of  the 


quadrilaterals.  If  the  mesh  has  also  triangular  elements  we  pro 
ceed  very  analogously. 


In  Step  1  we  use  the  mapping  M.  .  which  is  extended  on  S 

t  J 


and  consider  U.  as  the  image  of  u  on  G.  it  is  easy  to 

i  <  J  *  >  J 


show  that  the  extension  function  V(t,rf)  having  the  same  proper 
ties  as  mentioned  in  (6.7)  exists  for  T.  See  e.g.  [6]. 

All  other  steps  are  now  the  same  only  rendering  that  the 


"correction"  fucntions  now  could  be  of  degree  2p^  because 


is  polynomial  of  degree  2p  on  the  diagonal  of  S. 


7 .  Numer ical  Examples 

Let  us  consider  the  problem 

(7.1)  A  u  =  0  on  Q 

when  Q  is  an  L-shaped  domain  as  shown  in  Figure  7 . 1  and  the 
Dirichlet  conditions  are  prescribed  on  one  part  of  ^0  and  the 
Neumann  conditions  are  the  other  part  of  rfO . 


Figure  7.1.  The  domain  0. 

We  will  consider  two  problems  with  various  combinations  of  the 
Dirichlet  and  Neumann  boundary  conditions  with  the  exact  solution 
Case  A: 

1/3  1 

(7.1)  u  =  r  '  sin  ^ 

Case  B : 

2/3  2 

(7.2)  U  =  V  COS 

The  sequence  of  meshes  (<i  =  0.15)  is  characterized  by  the 


parameter  n  =  1 , . . . ,6  and  is  shown  in  Figure  7.2. 

Let  us  first  consider  the  case  where  the  Dirichlet  boundary 

condition  is  prescribed  on  the  entire  do .  Then  in  the  case  A  we 

have  g^  =  g  |  analytic  on  f.,  i  =  1,2,...  ,5  while  on 
i 

1/3 

r 6  :  gg  =  r  .  Now  we  can  use  Lemma  4.17  and  conclude  that  for 

2 

l  =  1,...,5,  g.  e  ©  ( f.),  0  <  ©  <  1  arbitrary  and  for  i  =  6  we 


.0225-7}  f.0225 


DETAIL 


0.15  =.0225 


Figure  7.2.  The  meshes  ft 


1 

6  ' 


Using 


obviously  have  i ■  =  ^  and  hence  g1  ^  ( r  6  )  with 

3/2  2 

Theorem  4.3  we  see  that  g  ^  B  (f)Q  )  with  <  >  ^ 

2  2  - 

Theorem  3.2  we  have  u  ^  ft  (Q)  with  3  >  maxl^,^ 


and  hence  by 

and  3  .  >  3  . 

1  1 


for  i  =  2,..., 6,  depend  on  the  problem.  In  our  case  it  can 

be  shown  that  3  =  -|  and  3.  =  0,  i  =  2,...,  5.  Hence  1  >  ■<  > 

1  v  1  ^ 

2 

—  arbitrary  and  0  <  3  .  <  1  arbitrary  for  i  =  2,..., 6.  This  of 
course  is  obvious  also  from  the  fact  that  the  solution  is  given  by 
(7.1).  Analogously  in  the  case  B  we  have  1  >  3  >  ,  and  0  < 


3  <  1  for  1  =  2 . 6. 

Denote  E(u)  =  -|b(u,u),  Epg(n,p)  =  E(Ug)  where  S  being 
characterized  by  the  mesh  Q™  and  degree  p  and 


B  ( u  ,  v )  = 


fdudv  <9u  dvl  ,  , 

1.53  55  +  5V  5^JlSxd'-' 


2 

Let  e  =  u-us<  i|ei!  =  E(e)  and  i|  e  !|  ^  =  n  e !!  £/ ij  u1!  £  be  the  error, 
the  error  norm  and  the  relative  error  norm,  respectively. 

We  have  in  the  case  A:  E(u)  =  0.423569  and  in  the  case  B: 

E ( u)  =  0.918113. 

Figure  7.3a,b  depict  the  relative  error  of  the  finite  element 
solution  in  0^  in  the  double  logarithmic  scale.  We  mention  that 
for  n  =  1  and  p  =  1  we  have  N  =  0  because  the  finite  element 
solution  is  determined  directly  by  its  values  at  the  boundary.  We 
see  that  in  the  case  A  the  p-version  does  not  practically  converge 
while  in  the  case  B  it  does.  Nevertheless  the  h-p  version 
(n  =  p)  converges  in  both  cases  as  the  theorems  5.1  predicts.  We 
also  show  in  the  figure  the  degree  of  the  elements.  Figure  7.4a,b 
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□ - □  n  =  5 
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Figure  7.3b.  The  relative  error  of  the  p  and  h-p  version  in 

i'e'i  „  '  Iq  N  scale.  Case  B. 


Figure  7.4a.  The  relative  error  of  the 

,1/3 


'N' 


p  and  h-p  version  in 
scale,  Case  A. 
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RELATIVE  ERROR  IN  THE  ENERGY  NORM  Hell 


CASE  B 


O.OI 


0.001 


X - X 

n  =  3 

\v 

A - A 

n  =4 

□ - □ 

n  =  5 

o — o 
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show  the  same  results  but  in  the  scale  €gi‘ eil_  *  N 
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that  the  error  of  the  h-p  version  is  nearly  linear  in  this  scale 


the  way  how  g  is  replaced  by  g  (see  Lemma  6.2). 
in  [9]  that  the  p-version  converges  with  optimal  rate  provided  that 
g  e  H1  (T  )  .  If  g  <z  H^(f)  to  get  optimal  rate  of  convergence  we 
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have  to  replace  g  by  g^  in  another  way  (so  called  H 
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jection,  see  [10]).  Then  the  convergence  is  also  guaranteed. 
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Figure  7.5  compares  the  performance  of  the  method  with  H  pro¬ 

jection  in  the  case  A  for  n  =  1.  We  show  in  Figure  7.5  the  slope 
-1/3 

°f  0(N  )  which  is  also  the  optimal  rate.  For  the  detailed 

comparison  of  the  performance  of  various  projections,  specifically 
1  1/2 

the  H  and  H  -projection  we  refer  to  [7]. 
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Figure  7.6a.  The  relative  error  of  the  h-p  version  for  various 
boundary  conditions,  Case  A. 
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In  Figures  7.6a,b  we  show  the  performance  of  the  h-p 
version  (p  =  n)  for  various  combinations  of  the  Dirichlet  and 
Neumann  conditions  (using  H* -pro jection  technigue  for  Dirichlet 
conditions).  We  see  that  there  is  no  significant  difference 
between  the  performance  of  the  h-p  version  for  various  combina¬ 
tions  of  the  boundary  conditions.  (We  mention  that  in  the  case  A 
the  Dirichlet  condition  on  T  ^  is  homogeneous  and  so  it  does  not 
contribute  to  the  error.)  In  contrast  the  p-version  with  the 
H1)!  )  projection  performs  independently  of  the  boundary  condi¬ 
tions  only  if  g^  H 1 ( r  ^ )  while  for  the  Dirichlet  condition  with 
-£  H° ( T  ^ ) ,  o  <  1  the  performance  deteriorates.  This  can  be  seen 
by  comparison  of  Figures  7.3a,b  resp.  7.7a,b  where  the  error  is 
given  for  the  Dirichlet  resp.  Neumann  boundary  condition.  We  see 
that  in  the  case  B  the  performance  of  the  p-version  (with  H1 
projection)  for  Dirichlet  boundary  conditions  is  the  same  as  for 
the  Neumann  condition  while  in  the  case  A  we  see  significant 
differences . 

If  the  Neumann  condition  or  Neumann  condition  and  homogeneous 
Dirichlet  conditions  is  prescribed,  then  the  strain  energy  of  the 
finite  element  solution  is  increasing  with  p,  i.e.,  E  (n,p  ) 
Epg(n,p2)  for  p2  >  p^.  Because  increasing  n,  the  shape  of 
elements  is  changed,  we  do  not  have  necessarily  E  (n  ,pi 
EFE<n2'P)  f°r  n2  >  ni  although  practically  this  usually  occurs. 
If  the  continuous  Dirichlet  condition  is  prescribed  on  the  entire 

boundary  and  g^  are  polynomials  of  degree  p  ?  pQ  then 
EpE(n,p2)  ?  Efn.p^  for  p2  >  px  -  p  ,  i.e.,  the  strain  energy  is 

decreasing  with  p.  If  the  Dirichlet  condition  is  not  a  polynomial 
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Figure  7.8  shows  the  behavior  of  E__  for  the  case  A,  where 

r  E 

n  =  1,2  and  p  =  1 . 8  when  the  monotonicity  occurs  only  in 

the  case  d  as  expected. 
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Figure  7.8.  The  monotonicity  of  the  energy  of  the  finite 
element  solutions. 


In  the  case  when  the  approximation  of  the  nonhomogeneous 
Dirichlet  boundary  conditions  do  not  contribute  to  the  error  we 


eg  =  |E(us)  -  E ( u ) 


In  the  case  when  the  Dirichlet  conditions  are  prescribed  on  the 
entire  dQ  we  have 

■e-!g  =  I  (E(us)  -  E(u)  )  +  *  I 


where  the  correction  term  k  is  due  to  the  fact  that  the  finite 
element  solution  does  not  satisfy  exactly  the  boundary  condition 


U  .  e  .  , 
and  I 


gg*g).  Nevertheless  k  is  usually  negligible  for  p  2 


E(us)  -  E(u) | 1/2  =  e 


is  very  close  to  e 


The  term  k 


can  be  easily  computed  if  the  small  solution  is  known.  In  fact 


e,  2  =  :u-us:  2  =  iB(us-u,us-u; 


B ( Us  » Ug )  +  B(U,U)  -  2B(u,us)j 


B ( us  * us )  -  B ( u , u ) 


+  B ( u , u-uc 


=  E  (  u. 


E  (  u  )  +  k 


where 


k  =  B(u,u-us) . 


Because  Au  =  0  we  have  by  integrating  by  parts 


•k  = 


C*U  /  i  j 

,7n  (U’US>  ds 


'dO 


and  u-u_  is  known  on  do  as  well  as  rr— .  Table  7.1a,b  shows 
b  a  n 

the  correction  term  k  for  the  mesh  0*},  n  =  6  and  the  relative 

o 


value  of  V/  e  .  depending  on  p  for  the  case  A  and  B. 


Table  7.1a.  The  correction  term  R 


CASE  A 


p 

R 

!le!|? 

E 

R/' !  e i!  2  9 

/  E  /o 

1 

-1 . 237 ( -2 ) 

5 . 626 ( -2  ) 

21.99(0) 

2 

-1 . 755 ( -4 ) 

4 . 945 ( -3 ) 

3.55  { 0 ) 

3 

-4 . 229 ( -6 ) 

3 . 419 ( -3 ) 

1.21(-1) 

4 

-1  .  203 ( -7 ) 

6 . 091 ( -4 ) 

1 . 97 ( -2 ) 

5 

1 . 053 ( -1 1 ) 

5 . 027 ( -4 ) 

2 . 09 ( -6 ) 

6 

3 . 927 ( -9 ) 

4 . 701 ( -4 ) 

8 . 35 ( -4 ) 

7 

4 . 064 ( -9 ) 

4 . 535 ( -4 ) 

8 . 96 ( -4 ) 

8 

4 . 059 ( -9 ) 

4 .433 ( -4 ) 

9 . 1 8 ( -4  ) 

Table  7.3b.  The  correction  term  R 
CASE  B 


P 

R 

'!  e.|  l 

E 

R/We\\l  J 

E  ° 

1 

-2 . 698 ( -2 ) 

5 . 08  2 ( -2 ) 

19.67(0) 

2 

-1 . 542 ( -4 ) 

1 . 695 ( -3 ) 

9.09  (0) 

3 

-2 . 614 ( -6) 

1  .  801 ( -4 ) 

1.44  ( 0 ) 

4 

-6. 1 26 ( -8 ) 

2  .  296 ( -5 ) 

2 . 67 ( - 1 ) 

5 

- 1 . 629 ( -9 ) 

3 . 336 ( -6 ) 

5 .04 ( -2  ) 

6 

-4 . 233 ( - 1 1  ) 

5 . 512 ( -7 ) 

7 . 68 ( -3 ) 

7 

7 . 636 ( - 1 2 ) 

1  .  202 ( -7 ) 

6 . 35 ( -3 ) 

8 

9. 2  68 ( - 1 2 ) 

4 . 6 1 9 ( -8  ) 

2 . 01 ( -3 ) 

Let  us  mention  that  for  higher  p  the  correction  is  influ 
enced  by  round  off  errors. 
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